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AP&TRACT 

A  microprocessor-controlled  magnetometer  which  accurately  locates  and  identifies 
compact  ferrous  objects  in  real-time  is  described.  The  person-portable  instrument  consists 
of  a  cart-mounted  cesium  magnetometer,  optical  encoder,  microcontroller,  interface  and 
laptop  computer.  The  instrument  guides  the  operator  to  collect  simultaneous  magnetic  field 
and  position  data  in  a  horizontal  plane  above  an  object.  Custom  algorithms  estimate 
location  and  dipole  moment  and  use  the  latter  to  classify  the  object.  Data  collection  takes 
6  -  13  minutes,  location  and  moment  estimation  5  seconds,  classification  30  seconds. 
Experiments  using  two  ferrous  spheroids  and  studies  using  magnetic  total  field  and  vertical 
component  magnetic  maps  generated  by  a  mathematical  computer  model  are  described. 
Limits  of  error  in  estimation  of  location  and  dipole  moment,  error  in  classification,  and 
relative  effects  of  sources  of  error  are  quantified.  The  rms  error  for  location  vector 
components  was  0.019  m  -  0.04f  rn  compared  to  the  average  precision  of  0.003  m  -  0.005  m. 
The  average  magnitude  of  the  difference  between  estimated  and  theoretical  dipole  moment 
vectors  as  a  percentage  of  the  theoretical  dipole  moment  was  24.5  ±11.4%  compared  to 
precision  uf  0.51  -  8.21%.  Pattern  classification  with  a  computer  generated  dipole  moment 
design  set  is  .  escribed.  Deviation  between  experimental  moment  estimates  and  the 
computer  model  degraded  performance  but  the  misclassification  rate  for  the  two  objects  for 
which  experimental  measurements  were  made  was  11.1%.  If  an  experimental  design  set 
were  used,  analysis  shows  that  the  limiting  misclassification  error  for  the  present 
experimental  precision  should  be  between  2  and  5%. 
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RfeSUMfe 

On  d6crit  un  magn^tomitre,  pilots  par  un  micropiocesseur,  qui  permet  de  localiser 
et  d’identifier  rapidement  en  temps  r6el  des  objets  ferreux  compacts.  L’appareil,  qui  pent 
etre  port^  par  une  personne,  est  constitu6  d’un  magn^tcm^tre  au  c6sium,  d’un  codeur 
optique,  d’un  micro  contrbleur,  d’une  interface  et  d’un  ordinateur  portatif,  le  tout  mont6  sur 
un  chariot.  L’appareil  permet  k  I’utilisateur  d’obtenir  simultan6ment  des  donn^es  sur  le 
champ  magndtique  et  la  position  dans  un  plan  horizontal  au-dessus  d’un  objet.  Des 
algorithmes  61abor6s  sur  mesure  permettent  d’estimer  le  moment  dipolaiie  et  sa  position 
et  d’utiliser  ces  donn6es  pour  classer  I’objet.  II  faut  de  6  it  13  minutes  pour  obtenir  les 
donndes,  5  secondes  pour  estimer  le  moment  et  ia  position  et  30  secondes  pour  classer 
I’objet.  On  d^crit  des  experiences  avec  deux  sph6roides  ferreux,  ainsi  que  des  6tudes 
effectu6es  avec  des  cartes  du  champ  magn6tique  total  et  de  la  composante  magn6tique 
verticale,  produites  ^  I’aide  d’un  module  mathematique  sur  ordinateur.  On  determine 
I’importance  des  limites  d’erreur  sur  I’estimation  du  moment  dipolaire  et  de  sa  position,  de 
I’erreur  de  classement,  ainsi  que  des  effets  relatifs  des  sources  d’erreur.  L’erreur 
quadratique  moyenne  sur  les  composantes  du  vecteur  de  position  etait  de  0,019  m  - 
0,045  m,  en  comparaison  avec  la  precision  moyenne  de  0,003  m  -  0,005  rn.  La  grandeur 
moyenne  de  I’ecart  entre  la  valeur  estimee  et  la  valeur  theorique  des  vecteurs  du  moment 
dipolaire,  exprimee  en  pourcentage  du  moment  dipolaire  th6orique,  etait  de  24,5  ±  1 1,4  %, 
en  comparaison  avec  la  precision  qui  etait  de  0,51  -  8,21  %.  On  decrit  le  classement  des 
resultats  realises  k  I’aide  d’un  ensemble  de  reference  de  moments  dipolaires  cree  par 
ordinateur.  L’ecart  entre  les  moments  estimes  experimentalement  et  les  valeurs  obtenues 
par  modeiisation  informatique  diminuait  le  rendement  de  I’appareil;  neanmoins,  le  taux  de 
classement  incorrect  pour  les  deux  spheroYdes  ayant  fait  I’objet  des  mesures  experimentales 
etait  de  11,1  %.  L’ analyse  revdle  que,  pour  un  ensemble  de  reference  experimental,  la 
limite  d’erreur  de  classement  devrait  Stre  comprise  entre  2  et  5  %  pour  la  pr6sente 
precision  expcrimentale. 
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Executive  Summary 

Magnetometers  have  been  used  to  detect  ferrous  objects  in  a  wide  range  of  military  and 
civilian  applications  since  the  1930’s.  Magnetometers  produce  a  signal  which  is  indicative 
of  magnetic  field  strength,  but  do  not  explicitly  provide  the  accurate  location  or  identity  of 
a  detected  object.  In  a  number  of  applications,  particularly  the  clearance  of  old  artillery 
ranges  and  mine  detection,  such  an  ability  would  be  very  desirable.  Several  research  groups 
in  several  countries  have  developed  data  logging  magnetometers  with  limited  location  and 
identification  capability  for  large  scale  surveys  of  magnetic  anomalies.  However,  all  the 
instruments  must  transfer  data  to  a  computer  which  analyses  the  data  off-line.  They  can 
localize  a  compact  ferrous  object  to  within  a  few  times  its  depth  of  burial  and  they  can 
classify  an  object  in  terms  of  its  rough  size.  None  of  them  can  explicitly  determine  the 
location  to  within  a  fraction  of  the  depth  and  explicitly  determine  the  identity  of  the  detected 
object  in  real-time. 

This  report  presents  an  improved  “smart”  magnetometer.  Based  on  an  earlier  prototype 
produced  by  this  laboratory,  the  new  version  is  the  first  magnetometer  having  the  ability 
to  explicitly  and  accurately  locate  and  identify  compact  ferrous  objects  in  real-time.  The 
instrument,  which  is  person-portable  but  could  be  modified  for  vehicle  mounting,  consists 
of  a  cesium  vapour  magnetometer  mounted  on  a  cart  with  a  wheel-mounted  optical  encoder, 
a  microcontroller,  interface  and  a  laptop  computer.  The  instrument  guides  the  operator  in 
the  collection  of  simultaneous  magnetic  field  and  position  data  in  a  horizontal  plane  above 
an  object.  Location  is  estimated  by  applying  a  custom  location  algorithm  to  the  data  and 
identity  is  established  by  a  custom  pattern  classification  using  a  byproduct  of  the  location 
algorithm,  the  dipole  moment,  as  a  feature  vector,  'fhe  instrument  is  more  robust  and 
user  friendly  than  the  earlier  prototype  and  its  accuracy  in  estimating  location  and  dipole 
moment  vectors  has  been  improved.  Design  sets  of  dipole  moments  (used  for  comparison 
in  the  pattern  classifier)  can  now  be  updated  in  the  field.  The  user  interface  uses  a  lap  top 
computer  to  communicate  with  the.  magnetometer’s  microcontroller  and  is  much  less  cryptic 
than  the  previous  version’s  keypad  and  four  digit  display.  Magi  tic  field  data  and  design 
sets  can  be  transferred  between  the  computer  and  the  rnicroc>  troller.  With  the  present 
instrument,  as  few  as  6  minutes  are  required  to  collect  the  data,  location  and  dipole  moment 
estimation  requires  5  seconds  and  classification  requires  30  seconds. 

To  quantify  the  penomiance  of  the  instrument,  a  detailed  study  was  recently  completed 
to  determine  the  error  in  estimating  the  location  and  dipole  moment  of  ferrous  spheroids. 
Spheroids  were  chosen  because  they  can  be  similar  in  shape  and  size  to  unexploded  ordnance 
and  have  magnetic  fields  that  well  approximate  those  of  unexploded  ordnance.  Also,  a 
mathematical  model  for  the  static  magnetic  field  of  a  spheroid  induced  by  a  homogeneous 
magnetic  field  exists  which  can  be  used  both  to  generate  design  sets  and  to  aid  in  analysing 
estimation  errors,  Experiments  were  performed  using  two  different  ferrous  spheroids. 

The  precision  in  estimating  location  varied  from  0.12  to  1.35%  of  the  depth  for  a  fixed 
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object  and  orientation  at  a  given  depth.  There  was  a  slight  variation  in  location  estimate  with 
object  type  and/or  orientation  which  increased  the  location  uncertainty  to  1.67  to  1.90%  of 
the  depth.  The  precision  in  estimating  the  dipole  moment  varied  from  0.51  to  8.21  %  of  the 
dipole  moment  magnitude.  There  was  a  slight  increase  in  uncertainty  as  depth  decreased, 
but  the  precision  was  not  closely  correlated  with  peak  magnetic  field. 

All  three  components  of  the  location  estimate  had  a  depth  dependent  bias  which  was 
larger  than  can  be  accounted  for  by  the  precision  of  the  estimation.  Errors  in  the  vertical 
component  were  generally  larger  than  those  for  the  horizontal  components.  The  root 
mean  square  (RMS)  error  over  all  depths  was  0.020  m  and  0.019  m  for  the  two  horizontal 
components  of  location  and  0.045  m  for  the  vertical  component.  By  comparison,  the 
average  precision  for  a  fixed  object  and  orientation  was  0.005  m  and  0.003  m  for  the  two 
horizontal  components  and  0.004  m  for  the  vertical  component. 

The.re  was  a  deviation  between  theoretical  values  and  experimental  estimates  of  the 
dipole  moment  that  was  greater  than  the  uncertainties  in  the  estimates.  The  unweighted 
average  difference  between  estimated  and  theoretical  dipole  moment  components  as  a 
percentage  of  the  theoretical  dipole  moment  was  4.8  ±  7.6%  and  -3.3  ±  15.5%  for  the  two 
horizontal  components  of  the  dipole  moment,  -10.8  ±  17.5%  for  the  vertical  component 
and  24.5  ±  1 1 .4%  for  the  magnitude  of  the  vector  difference.  There  was  no  clear  trend  with 
object  type,  orientation  or  depth. 

Following  quantification  of  the  error  in  the  estimation  of  location  and  dipole  moment, 
sources  of  error  and  their  relative  effects  were  analysed.  By  applying  the  location  algorithm 
to  magnetic  total  field  and  vertical  component  data  generated  by  mathematical  computer 
model,  the  relative  contribution  of  the  various  error  sources  could  be  estimated.  It  was 
argued  that  a  substantial  fraction  of  the  uncertainty  (precision)  in  dipole  location  and 
moment  vector  estimates  was  due  to  positional  uncertainty  in  the  experiments.  It  was  also 
demonstrated,  however,  that  positional  error  in  these  experiments  had  little  influence  on 
the  overall  location  and  moment  vector  estimation  and  that  remnant  magnetization  in  the 
spheroids  must  be  small  compared  to  the  induced  magnetization.  It  was  shown  that  the 
finite  volume  of  the  sensor  head  was  not  a  significant  source  of  error  for  the  geometry  of 
these  experiments.  A  variation  of  error  with  depth  and  obi'ct  orientation  was  noticed  for 
the  computer  generated  data.  This  suggested  that  the  main  contributor  to  the  error  was 
estimation  of  the  position  of  the  field  extrema  and  the  use  of  total  field  in  place  of  vertical 
Component  magnetic  data.  Ultimately,  it  was  estimated  that  the  error  due  to  the  total  field 
approximation  varied  from  14  to  18%  of  the  moment  magnitude  and  that  the  error  due 
to  estimation  of  extrema  positions  and  higher  order  moments  varied  from  0  to  8%  of  the 
moment  magnitude. 

The  analysis  of  error  in  the  computer  generated  total  magnetic  field  data  also  allowed 
limits  to  be  set  on  the  mi.sclassification  rate  that  can  ultimately  be  attained.  The  moment 
estimates  for  a  given  object  and  orientation  were  found  to  vary  by  between  2  and  4%  of 
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the  moment  magnitude  as  the  depth  varied.  If  the  estimated  moments  were  used  in  the 
design  set,  this  suggests  that  the  limiting  misclassificatiori  error  would  be  between  roughly 
1.4  and  2.4%  if  the  experimental  (mainly  positional)  error  v/ere  substantially  less  than  the 
algorithmic  and  approximation  errors.  For  the  present  experimental  precision  (4.1  ±  3.0% 
),  the  misclassification  eiror  should  be  roughly  between  2  and  5%  (assuming  quadrature 
error  summation)  if  the  design  set  were  experimentally  obtained. 

Pattern  classification  was  performed  with  a  computer  generate  dipole  moment  design  set 
consisting  of  8  objects  including  the  2  used  in  these  experiments.  The  gross  misclassification 
rate  was  10.3%  for  the  class  number  4  spheroid  and  100%  for  class  number  6.  Closer 
examination  revealed  that  class  6  was  always  classified  as  class  number  3,  which  is  a 
similarly  shaped  spheroid.  This  is  likely  due  to  the  deviation  between  the  moment  estimates 
for  the  experimental  spheroid  and  the  computer  model.  The  repeatability  of  classification, 
that  is,  the  percentage  of  cases  in  which  an  object  is  classified  as  the  same  class,  is  a  better 
measure  of  classifier  performance.  The  repeatability  of  classification  was  91.7%.  If  we 
consider  only  the  two  objects  for  which  experimental  measurements  were  made,  we  see 
that  only  4/36  cases  (11.1%)  were  classified  incorrectly. 

The  overall  performance  of  the  smart  magnetometer  is  very  encouraging.  The  report 
closes  with  a  discussion  of  future  work  that  must  be  done  to  improve  the  instrument  and  to 
make  it  successful  as  a  practical  locator  and  identifier  for  buried  ferrous  ordnance. 
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1.  Introduction 


Magnetometers  have  been  used  to  detect  ferrous  objects  in  a  wide  range  of  applications 
since  the  1930’s.  Magnetometers  produce  a  signal  which  is  indicative  of  magnetic  field 
strength,  but  do  not  explicitly  provide  the  accurate  location  or  identity  of  a  detected  object. 
In  a  number  of  applications,  particularly  the  clearance  of  old  artillery  ranges  and  mine 
detection,  such  an  ability  would  be  very  desirable.  Several  research  groups  have  addressed 
the  problem  of  large  scale  surveys  to  detect  magnetic  anomalies.  A  vehicle-towed  detector 
array,  STOLS,  has  been  developed  by  Geo-Centers  Inc.,  Newton  Falls,  MA,  USA,  for  the 
US  Navy  [1]  to  locate  magnetic  anomalies  over  areas  of  several  thousand  square  meters. 
The  device  consists  of  a  set  of  commercially  available  cesium  vapour  magnetometers  and 
a  microwave  tiiangulation  position  measurement  system  connected  to  a  data  collection 
system.  Data  is  processed  off-line  to  produce  field  intensi  ty  maps  which  can  be  interpreted 
to  roughly  locate  magnetic  objects.  Further,  compact  magnetic  objects  such  as  artillery 
shells  can  be  located  more  accurately  and  roughly  grouped  in  size,  based  on  dipole  strength, 
using  a  fairly  slow  iterative  nonlinear  least  squares  fit  dipole  locator  algorithm  originally 
developed  by  our  laboratory  [2 1 .  Another  large  area  magnetic  field  mapping  instrument,  the 
TM-3,  has  been  developed  by  the  University  of  New  England,  Australia  [3].  This  hand-held 
instrument,  which  is  based  on  a  single  cesium  vapour  sensor,  uses  a  data  logger  to  collect  a 
magnetic  field  map.  Data  is  transferred  to  a  PC  which  produces  images  of  the  magnetic  field 
in  a  plane  and  allows  rough  localization  and  identification  of  anomalies  by  visual  inspection 
of  the  shape  of  the  magnetic  field.  Institut  Dr.  Forster,  Reutlingen,  Germany,  manufactures  a 
hand-held  fiuxgaie  gradiometer,  the  FEREX  CAST  4.021.06,  which  is  intended  for  smaller 
areas.  It  can  log  data  while  scanning  an  area,  then  transfer  the  data  to  a  PC  for  off-line 
analysis  [4].  A  similar  data  logging  magnetometer  with  a  PC  analysis  package,  called 
CAMAD,  was  announced  some  time  ago  by  Aprotec  Ltd.,  Manchester,  UK  [5],  but  nothing 
has  been  heard  of  it  since  the  initial  announcement.  All  these  instruments  must  transfer 
data  to  a  computer  which  analyses  the  data  off-line.  They  can  localize  a  compact  ferrous 
object  to  within  a  few  times  its  depth  of  burial  (the  exception  being  STOLS  using  the  DRES 
iterative  location  algorithm)  and  can  classify  an  object  in  terms  of  its  rough  size.  None  of 
them  can  explicitly  determine  the  location  to  within  a  fraction  of  the  depth  and  explicitly 
determine  the  identity  of  the  detected  object  in  real-time. 

In  [6],  a  real-time  method  to  explicitly  estimate  the  location  and  determine  the  identity 
of  a  compact  ferrous  object  was  presented.  It  was  based  on  sampling  the  magnetic  field 
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of  the  object  in  a  horizontal  plane,  storing  simultaneous  position  and  magnetic  field  data, 
and  then  using  a  noniterative  algorithm  to  determine  the  location  and  components  of  the 
dipole  moment  associated  with  the  object.  A  byproduct  of  the  location  algorithm  was  an 
estimate  of  the  dipole  moment  which  was  used  by  a  DRES  continuous  parameter  (CP) 
pattern  classifier  to  identify  the  object.  A  prototype  microprocessor-controlled  cesium 
vapour  “smart”  magnetometer  was  described  which  collected  simultaneous  magnetic  and 
position  data  and  then  used  these  algorithms  to  locate  and  identify  ferrous  spheres  and 
spheroids.  A  few  preliminary  results  were  also  reported.  This  was  the  first  magnetometer 
with  the  ability  to  explicitly  locate  and  identify  compact  ferrous  objects  and  the  first  total 
field  magnetometer  with  explicit  object  location  estimation  in  a  self-contained  instrument. 

This  report  presents  an  improved  “smart”  magnetometer.  The  new  version  is  the  first 
magnetometer  having  the  ability  to  explicitly  and  accurately  locate  and  identify  compact 
ferrous  objects  in  real-time.  The  instrument,  which  is  person-portable,  consists  of  a  ce¬ 
sium  vapour  magnetometer  mounted  on  a  cart  with  a  wheel-mounted  optical  encoder,  a 
microcontroller,  inteiface  and  a  laptop  computer.  The  instrument  guides  the  operator  in 
the  collection  of  simultaneous  magnetic  field  and  position  data  in  a  horizontal  plane  above 
an  object.  Location  is  estimated  as  before  by  applying  the  DRES  noniterative  location 
algorithm  to  the  data  and  identity  is  established  by  pattern  classification  using  the  dipole 
moment  as  a  feature  vector.  'Fhe  instrument  is  now  more  robust  and  user  friendly  than 
the  earlier  prototype  and  its  accuracy  in  estimating  location  and  dipole  moment  vectors 
has  been  improved.  Design  sets  of  dipole  moments  (used  for  comparison  in  the  pattern 
classifier)  can  now  be  updated  in  the  field.  The  user  interface  uses  a  lap-top  computer  to 
commmunicaie  with  the  magnetometer’s  microcontroller  and  is  much  less  cryptic  than  the 
previous  version’s  keypad  and  four  digit  display.  Magnetic  field  data  and  design  sets  can 
be  transferred  between  the  computer  and  the.  microcontroller. 

To  quantify  the  peifonnance  of  the  instrument,  a  detailed  study  was  recently  completed 
to  determine  the  error  in  estimating  the  location  and  dipole  moment  of  ferrous  spheroids. 
Spheroids  were  chosen  because  they  can  be  similar  in  shape  and  size  to  unexploded  ordnance 
and  have  magnetic  fields  that  well  approximate  those  of  unexploded  ordnance  [7],  [8].  Also, 
a  mathematical  model  for  the  static  magnetic  field  of  a  spheroid  induced  by  a  homogeneous 
magnetic  field  exists  [8]  which  can  be  u.sed  both  to  generate  design  sets  and  to  aid  in 
analysing  estimation  errors.  The  present  study  has,  in  fact,  identified  several  sources  of 
error  in  an  attempt  to  leain  if  and  how  they  might  be  ameliorated. 

Chapter  2.  provides  the  necessary  theoretical  framework  regarding  magnetic  field  mea¬ 
surements,  the  locatio.n  and  moment  estimation  algorithm  and  the  pattern  classifier.  Also 
the  mathematical  model  for  the  magnetic  field  and  dipole  moment  induced  in  a  ferrous 
spheroid  by  a  uniform  magnetostatic  field  is  also  developed.  The  model  is  needed  to  gen¬ 
erate  design  .sets  for  the  classifier  as  well  as  to  analyse  the  performance  of  the  location 
algorithm.  Chapter  3  describes  the  magnetometer,  the  experimental  layout,  the  data  collec¬ 
tion  procedure  and  the  initial  calibration  of  the  instrument.  Chapter  4  presents  the  results 
of  applying  the  location  estimation  algorithm  to  a  large  number  of  magnetic  field  maps 
of  different  objects,  orientations  and  depths.  In  Chapter  5,  the  errors  associated  with  the 
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location  algorithm  are  analysed.  This  is  facilitated  by  comparison  with  location  estimation 
results  based  on  total  magnetic  field  maps  and  vertical  component  magnetic  firld  maps 
;  generated  by  computer  using  the  mathematical  model.  Performance  of  the  pattern  classifier 

is  discussed  in  Chapter  6.  Chapter  7  summarizes  the  performance  of  the  instrument  and 
discusses  possible  improvements  to  the  instrument  and  procedures  and  further  studies  that 
i  should  be  done. 

1  • 
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2.  Theory 

2.1  Magnetic  Field  Measurement 


The  spatial  extent  of  the  measureable  magnetic  fields  associated  with  most  compact 
objects  of  interest  is  generally  less  than  a  few  hundred  meters  and  the  time  taken  to  measure 
such  fields  is  less  than  a  few  minutes.  Fortunately,  the  earth’s  field  is  constant  over  such 
distances  (magnetic  induction  gradient  ~  10  nT/km)  and  does  not  change  significantly 
during  this  time.  In  spite  of  occasional  field  nonuniformities  arising  from  local  magnetic 
phenomena,  one  can  in  practice  usually  assume  that  the  ambient  field  is  constant.  A  detailed 
description  of  the  spatial  and  temporal  variation  of  the  earth’s  magnetic  field  and  the  field 
due  to  common  terrestrial  sources  may  be  found  in  [9]. 

ITierc  are  two  types  of  magnetometers  -  vector  sensors  and  total  field  sensors.  The  latter 
includes  self-oscillating  optically  pumped  cesium  vapour  magnetometers,  such  as  the  one 
used  in  the  present  work.  Total  field  magnetometers  measure  the  magnitude  of  the  field  but 
not  its  direction.  Their  chief  advantage  over  vector  sensors  is  that  the  fonner  arc  insensitive 
to  small  changes  in  the  orientation  of  the  sensor.  Principles  of  operation  of  cesium  vapour 
magnetometers  and  low  noise  measurement  techniques  may  be  found  in  [9]. 

We  will  now  discuss  the  signal  measured  by  a  total  field  magnetometer.  ITie  sec¬ 
ondary  field  (field  due  to  the  object)  in  a  cartesian  coordinate  system  is  denoted  by 
B  =  {B\,  Bz,  Bz)^',  where  the  superscript  T  denotes  the  transpose.  T the  primary  (ambient, 
i.e.,  earth’s)  field  is  Bo  =  {Bo\ ,  i?o2,  then  a  total  field  magnetometer  would  measure 


B,^- 


Bo  +  /5||  =  +  B^  +  2Bo  •  b)^''^  . 


Generally  Bo  >  B  and  a  Taylor  expansion  can  be  perfonned  to  obtain 


(2.1) 


.Bj  i?0  +  7  •  ^  + 


2Bo 


B^ 


where 

l^Bo^Bo.  (2.3) 

To  first  order,  a  total  field  magnetometer  measures  the  magnitude  of  the  earth’s  field  plus 
the  projection  of  the  magnitude  of  the  object’s  field  along  the  direction  of  the  earth’s  field 
vector. 
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2.2  Estimation  of  Dipole  Location  and  Moment 

A  space-fixed  cartesian  coordinate  system  with  arbitrary  origin  is  assumed  (Fig.  2.1). 
Assume  also  a  cartesian  coordinate  system  centered  at  the  dipole  location  (called  the  “dipole 
system”).  Vector  and  tensor  components  in  the  space-fixed  coordinate  system  arc  denoted 
by  upper  case  symbols  and  those  in  the  dipole  system  arc  denoted  by  lower  case  symbols. 
The  .system  is  oriented  such  that  m  -  {m\ , 0,  i.c.,  m2  =  0  and  the  X3  axis  is  parallel 

with  the  .V3  axis.  The  X2  axis  makes  an  angle  a  with  the  X2  axis.  Field  measurements  are 


Geometry  for  magnetostatic  dipole  location. 


made  in  the  X3  =  xo3  plane.  The  measurement  plane  is  chosen  to  be  at  X3  =  0,  without  loss 
of  generality,  and  A03  =  -X03  is  referred  to  as  the  “depth”  of  the  dipole.  The  measurement 
plane  is  normally  horizontal  and  the  X3  axis  is  normally  vertical.  This  is  the  geometry  used 
in  the  experiments  of  this  report  and  is  similar  to  that  found  in  most  practical  unexploded 
ordnance  detection  applications. 

Various  analytical  and  numerical  methods  have  been  derived  to  estimate  the  six  inde¬ 
pendent  dipole  parameters  (the  three  components  of  the  location,  Xo,  and  moment,  M) 
from  measurements  of  the  dipole  field  or  gradient  [61.  The  method  used  in  this  study  is  an 
algorithm  developed  at  this  laboratory  which  we  have  dubbed  the  extremum  method.  It  has 
the  advantage  that  it  is  noniterative,  which  makes  it  robust  and  quick.  Most  imix)rtaritly,  it 
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is  the  only  noniterative  algorithm  for  a  generally  oriented  dipole  which  can  use  input  data 
from  a  total  field  magnetometer.  The  algorithm  is  outlined  below  and  is  described  in  detail 
in  [10]. 

If  both  a  maximum  and  a  minimum  of  a  field  quantity  can  be  found  in  the  plane  of 
measurement,  six  independent  pieces  of  information  are  available  and  the  dipole  parameters 
can  in  principle  be  determined.  Bi  has  more  than  one  extremum  and  hence  is  a  candidate 
for  such  a  method.  is  not  itself  a  rotationally  invariant  quantity.  However,  if  the  angle 
between  the  normal  to  the  measurement  plane  and  earth’s  field  direction  is  <~  1 5  -20"  (e.g., 
high  to  mid-latitudes),  the  3-component  of  the  secondary  field  is  reasonably  approximated 
by  (i?t  -  Bq)  (Equation  2.2)  and  hence  B^  can  be  measured  by  a  total  field  magnetometer 
which  is  insensitive  to  orientation.  Note  mat  the  angle  in  question  at  our  laboratory  site  is 
17". 


It  can  be  shown  [  10]  that  the  positions  (rtf  ^ ,  xl- )  of  the  extrema  of  B^  in  the  measurement 
plane  are  given  by 

=  ;  rtf,- -0  i  =  1,2,3  (2.4) 

where  are  monotonic  analytic  functions  of  /?  =  m\jmz  whose  ranges  are  2  <  /i  <  oo, 
0  <  /a  <  1/2  and  —2<h<  —1/2.  The  positions  of  the  extrema  lie  in  the  plane  along 
a  straight  line  which  coincides  with  the  rta  axis.  (An  exception  is  the  case  of  a  vertical 
dipole  for  which  choice  of  rta  cxis  is  arbitrary.  A  maximum  exists  immediately  above  the 
dipole  and  a  continuous  ring  of  minima  is  centered  on  the  maximum.)  Thus,  the  angle  a 
which  relates  the  orientation  of  the  dipole  coordinate  system  to  the  space-fixed  system  is 
immediately  known  once  the  positions  of  the  extrema  arc  found.  It  turns  out  that ;  =  2 
corresponds  to  a  field  maximum  while  ;  =  3  is  a  field  minimum  B^'^^  and  j  =  1  is 
a  saddle  point.  The  field  values  at  the  maximum  or  min -mum  are  given  by 


n 

B,  =  — m;ro3 


o/  (l  +  (l  -h  [-1  +  3  (1  +  /J)"'  (1  +  /3/;) 


(2.5) 


where  m  =  and  /  =  2  (maximum)  or  3  (minimum). 

The  ratio  of  the  field  minimum  to  the  field  maximum  is  a  monotonic  function  of  P  which 
can  be  inverted  by  approximating  the  function  for  intervals  of  p  by  polynomials  of  the  form 

^  =  (2-6) 
fc=0 


where  ajk  are  the  fitted  coefficients  for  the  jth  interval  of  p.  Three  intervals  with  3  < 
n  <  6  for  each  interval  have  given  very  accurate  approximations  to  the  function.  The 
corresponding  coefficients  arc  given  in  [10]. 

The  algorithm,  then,  is  simple  to  implement.  First  the  space-fixed  coordinates  of  the 
position  in  the  measurement  plane  of  the  maximum  and  minimum  (Xf^,  X^j)  are  found 
and  a  is  estimated.  The  differences  in  .ti  coordinates  between  the  maximum  and  minimum 
are  obtained  by  a  simple  rotation  transformation  using  n.  Next,  P  is  deduced  from  the 
measured  ratio  of  field  minimum  to  maximum  by  Eq.  2.6  and  is  used  to  estimate  /z,  h- 
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Estimates  of  /2,  /a  and  the  differences  in  coordinates  of  maximum  and  minimum  are 
inserted  in  Eq.2.4  to  estimate  AToa.  The  xi  coordinates  of  the  maximum  and  minimum  are 
next  calculated  using  Eq.2.4  and  are  used  to  find  A'oi ,  X02  from, 

A'oi  =  A'fj  -  x\j  cos  a  ;  AToj  =  sin  a  (2.7) 

where  j  =  2, 3.  The  position  ATo  is  now  known  and  since  P,  a  are  known,  the  direction  of 
M  is  completely  specified.  Finally,  m,  the  magnitude  of  A/,  is  obtained  from  Eq.2.5. 

Small  modifications  to  the  algorithm  as  desciibed  are  necessary  for  the  case  of  a  vertical 
dipole,  and  some  heuristics  are  used  to  handle  minima  that  are  too  small  to  be  observed 
[10]. 


2.3  Identification  of  a  Spheroid  From  Its  Dipole  Moment 


The  previous  algorithm  can  estimate  the  locatiM  and  moment  of  the  dipole  associated 
with  a  compact  ferrous  object.  However,  it  is  the  location  and  identity  of  the  compact 
ferrous  object  that  is  required.  Six  parameters  uniquely  define  dipole  source,  but  with 
the  exception  of  the  sphere,  more  than  six  parameters  are  needed  to  specify  the  compact 
object.  This  means  that  if  field  measurements  are  made  at  suffxient  distance  from  a  source 
so  that  higher  order  multipoles  are  negligible  with  respect  to  ihe  dip<ole,  one  cannot  in 
general  distinguish  between  two  compact  orientable  bodies.  Information  is  necessary  from 
higher  order  multipoles  if  the  problem  is  to  be  have  a  unique  inverse.  Estimation  of  these 
multipoles  is  difficult,  usually  relying  on  fitting  the  field  to  a  source  model  and  solution  is 
dependent  on  the  shape  of  object. 

Fortunately  in  practical  applications  there  are  usually  a  small  number  of  object  shapes 
and  si/xs  applicable  to  a  particular  problem.  If  the  dipole  field  associated  with  each  object 
of  the  set  is  sufficiently  different  from  that  of  the  other  objects,  identity  can  be  reliably 
determined  from  the  dipole  field.  We  will  now  outline  the  method  that  is  used  by  the  smart 
magnetometer  to  identify  compact  ferrous  objects  based  on  their  magnetic  dipole  moments. 
It  employs  a  novel  pattern  classification  technique  which  is  discussed  more  fully  in  [11]. 


It  is  assumed  that  we  have  determined  the  location  and  dipole  moment  components  of 
the  object  relative  to  a  space-fixed  ca  ’csian  coordlriatc  sysietii  by  the  method  described 
in  Section  2.2  and  that  the  dipole  location  is  at  or  near  the  geometric  center  of  the  object. 
The  latter  turns  out  to  be  a  reasonable  assumption  (see  discussion  of  the  sources  of  dipole 
parameter  estimation  errors  in  Section  5.1).  ITie  plan  is  to  use  a  pattern  classification 
approach  by  which  an  object  under  scrutiny  is  identified  by  comparing  a  vector  composed 
of  features  (“test  vector”)  with  feature  vectors  from  a  set  of  known  objects  (“design  set”). 
The  feature  vector  chosen  will  be  the  dipole  moment  vector. 


The  space-fixed  components  of  the  dipole  moment  induced  in  a  homogeneous,  perme¬ 
able,  axially  symmetric,  compact  object  by  a  unifonn,  static  external  magnetic  field  will 
vary  with  the  orientation  of  the  object.  Since  two  angles  define  the  orientation,  the  locus 
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of  all  possible  dipole  moments  for  such  an  object  is  actually  a  two  dimensional  surface  in 
a  three  dimensional  space.  (The  dipole  moment  is  usually  independent  of  the  magnetic 
permeability  of  a  ferrous  object  in  the  earth’s  held  and  so  only  one  surface  exists  for  each 
unique  object  shape  and  size.)  For  a  spheroid,  which  also  has  .^ore-aft  symmetry,  unique 
values  of  the  space-fixed  magnetic  dipole  moment  M  occur  only  for  the  following  range 
of  the  two  continuous  orientation  parameters  (9,  <l>)  (6  =  (j>  ~  0),  {0  <  6  <  tt/I  when 
0  <  <t>  <  Itt)  and  {6  —  7r/2  when  0  <  <^  <  jt). 

Assume  that  a  design  set  of  dipole  moments,  Mj^k,  have  been  nt(.asiue?l  at  ilisL;cete 
intervals  (Oj,<j)k)  over  the  range  of  angles  for  which  unique  values  of  M  occur.  Consider 
the  region  of  the  dipole  moment  surface  corresponding  to  class  i  for  which  9j  <  0  <  ^j+i 
andi^jt  <  (f>  <  4>k+\-  This  region  of  the  surface  may  be  approximated  by  two  triangles.  One 
triangle  passes  through  points  Mj,k',  Mj+\  y,  Mj^k+i  (index  i  is  suppressed)  and  is  bomided 
by  Uj,k,  Vj^k  and  Uj,*..  -  Vj,k,  where 

Uj,k  =  -  ^j.k  ,  Vj,k  =  Mj+i,k  —  ^j.k‘  (2.8) 

The  other  triangle  passes  through  points  Mj+iy,  Mj+\,k+i ;  Mj,k+i  and  is  bounded  by  u '  j.- 
where 

If  X  is  a  test  vector  (suppressing  subscripts  j,  k  where  unambiguous  to  do  so),  then 

yj,k  =  ^-Mj.k-  (2.10) 

If  y  is  he  projection  of  y  onto  the  unprimed  mangle,  a  Gram-Schmidt  construction  gives 

=  Pi,.i,ku  +  qi,j,kv  (2.11) 

where  Pi,j,k,<ii,j.k  are  scalars  given  by 


(7.  ,■  i.  = 


(u'^'uj  -  [y'‘n)  {n  'v) 
{v  '^v)  (u  '^u)  —  (u'-'  v)^ 


Pi,j,k 


(y  '^uj  -  q 


1  0 


If  f  is  a  sample  from  the  class  i  corresponding  to  the  region  of  the  dipole  moment  surface 

"a  ^ 

bounded  by  v,vSl  -  v,  then  estimates,  0  and  <j),  of  the  continuous  parameters  associated 
wi^h  X  may  be  obtained  from 


0  =  Oi^qi,i,k{O,^x-9,)  ,  (2.13) 

1  =  <|>k^-  P,,j,k  i<Pk+i  -  <l>k)  ■  (2.14) 
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The  minimum  distance,  from  the  test  vector  to  the  triangle  is  given  by 

■  (2.15) 

The  previous  equations  give  the  minimum  distance  to  the  triangle  if  lies  within  the 
boundaries  of  the  triangle.  This  is  true  provided  p,,j,A,  >  0,9.,;,fe  >  OandO  <  Pi,j,k  +  qi,j,k  < 
1 .  If  these  conditions  are  not  satisfied,  then  is  replaced  by  the  minimum  distance  from 
the  test  vector  to  the  line  segments  which  bound  the  triangle. 

A  similar  process  is  applied  to  the  primed  triangle,  yielding  a  minimum  distance  d  ( 
Ihe  minimum  distance,  di,  from  the  test  vector  to  the  class  i  dipole  moment  surface  is  then 
approximated  by 


di  =  min{/i.j.fc,d'^}.  (2.16) 

The  dipole  moment  surface  for  d  =  0  is  independent  of  (^.  In  this  case,  the  manifold 
between  0i  ^  0,  0i+\ ,  <pk  and  (f>k+i  is  approximated  by  a  single  triangle  connecting  Mi, 
Mi^\,k 

This  classification  method  has  been  tested  on  a  design  set  consisting  of  noise-free 
magnetic  moments  for  six  different  spheroids,  typical  of  the  size  and  shape  of  a  wide 
variety  of  artillery  shells  (0.01  <  a  <  0.09  m,  2.5  <  e  <  3.75,  where  a,  c  are  defined  in 
Section  2.4).  These  were  computer  generated  at  15°  increments  of  the  orientation  angles 
(j).  Test  vectors  were  computer  generated  at  5°  increments  with  additive  Gaussian  noise 
of  different  levels.  The  probability  of  misclassification  was  about  1%  for  noiseless  test 
vectors,  3%  for  moments  with  noise  which  was  5%  of  the  moment  component  value  for 
each  component,  and  7%  for  moments  with  10%  noise  [11]. 

The  identification  method  has  been  presented  for  a  compact  axially  symmetric  object 
but  can  be  generalized  to  a  compact  body  of  arbitrary  shape  [11]. 


2.4  Multipole  Expansion  of  the  Magnetic  Field  of  a  Spheroid 

In  this  Section,  we  present  a  mathematical  model  for  the  magnetic  Reid  induced  in  a 
homogeneous  spheroid  of  arbitrary  orientation  by  a  uniform  magnetostatic  field.  The  model 
can  be  used  both  to  develop  a  design  set  for  use  by  the  magnetometer’s  pattern  classifier  (see 
previous  Section)  and  to  provide  controlled  magnetic  field  data  for  analysis  of  the  location 
estimation  algorithm  error.  The  model  has  been  shown  to  provide  magnetic  fields  that  are  a 
good  approximation  to  those  of  real  spheroids  and  some  unexploded  ordnance  in  the  earth’s 
field  [7],  [8].  A  brief  outline  of  the  model  is  given  below.  A  detailed  derivation  may  be 
found  in  [9]. 

It  is  assumed  that  a  homogeneous  uniformly  permeable  spheroid,  with  magnetic  perme¬ 
ability  Hr\Po  and  no  permanent  magnetization,  sii::  in  a  homogeneous  surrounding  medium 


Ur  ^CLASSIFIED 


DRES-SR-582 


with  permeability  /Xr2/<o.  where  fio  is  the  permeability  of  free  space  (Fig.  2.2).  There 


Figure  2.2 

Geometry  for  calculation  of  the  magnetic  field  of  a  spheroid. 

exists  a  uniform  parallel  external  magnetostatic  field.  A  body-fixed  cartesian  coordinate 
system  whose  3-axis  coincides  with  the  symmetry  axis  of  the  spheroid,  has  an  origin  at  the 
geometric  center  of  the  spheroid.  There  is  also  a  space-fixed  coordinate  system  to  which 
measurements  are  referenced.  Without  loss  of  generality,  the  space-fixed  system  is  chosen 
SO  that  the  primary  magnetic  field,  Bo,  has  no  2-component,  i.c.,  Bq  =  (J9oi ,  0,  Boa)  •  The 
spheroid  center  is  located  at  Xq  =  (A'^oi ,  A'oz,  Xoa)^  in  this  system.  Vectors  are  indicated  by 
lower  case  letters  in  the  body-fixed  system  and  upper  case  letters  in  the  space-fixed  system. 
The  angle  between  the  spheroid  symmetry  axis  and  the  X3  axis  is  0.  The  angle  between 
the  projection  of  the  symineU7  axis  on  the  horizontal  plane  and  the  Xi  axis  is  <l>.  It  can 
be  shown  [12]  that  the  induced  field,  and  hence  the  magnetization,  inside  the  spheroid  is 
uniform  and  parallel.  By  carrying  out  a  multipolc  expansion  of  the  static  field,  and  with  the 
assumptions  of  unifonn  parallel  internal  magnetization,  it  may  be  shown  that  the  secondary 
induced  magnetic  induction  measured  at  a  jjoint  r  =  (xi ,  X2,  xj)  in  free  space,  complete  to 
fourth  order,  is  given  by  |9] 

b^b^)  +  bW  (2.17) 

where  h^^'>  is  a  dipole  field  term  and  is  an  octupole  field  term.  The  two  field  terms  may 
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be  expressed  in  component  notation  using  summation  convention  as 

4®  =  +  3r-^  [x,m<‘>]  X.)  (2.18) 

and 

fea  ^  +  35r~'^x^xpx^xsrn}'^}^^ 

(2.19) 

where  is  the  dipole  moment  vector,  and  m^oli  the  rank  3  octupole  moment  tensor. 
Note  that  there  is  no  monopole  field  term,  as  expected,  and  that  the  quadrupole  term  is 
missing.  The  disappearance  of  the  latter  and,  in  fact  all  even  multipolc  moments 
n  an  integer),  derives  from  the  restrictions  on  the  magnetization  of  the  body.  A/',  together 
with  the  axial  and  fore-aft  symmetry  of  the  spheroid. 

The  dipole  moment  is  given  by 

=  a  =  1,2,3  (2.20) 

where  V  is  the  spheroid  volume.  For  the  octupole  moment  tensor,  only  10  of  the  27 
elements  are  independent  and  axial  symmetry  reduces  that  number  to  6.  These  are 

mffl  =  37ngl  =  3M;/„  ;  = 


(3)  (3) 

=  m)j3  = 

;  mgl  =  A/, 7,3  ; 

rnfi  =  M^Iys 

(2.21) 

where 

.  -  xf)  dv‘  ;  /33  =  j 

xfdF  . 

(2.22) 

Note  also  that 

rn-olp  = 

(2.23) 

"‘S7  =  ^  a  -7^  ^  7  ^  Of  .  (2.24) 


'Fhe  analysis  to  this  point  applies  to  any  axially  symmetric  body  with  fore-aft  symmetry. 
We  now  assume  that  the  body  is  a  spheroid  whose  symmetry  axis  length  is  2ae  and  whose 
maximum  diameter  orthogonal  to  the  symmetry  axis  is  la.  It  is  simple  to  show  that 


V'  =  -"Kt-O? 
3 


r  ^TT  5 

/„  = 


h,  =  -jje  u 


(2.25) 


Assume  that  the  external  magnetic  field  (usually  the  earth’s)  in  the  absence  of  the 
spheroid  is  60  =  (f>oi ,  kn-  ^03)^-  A  solution  of  the  boundary  value  problem  vields  [12] 

;•  =  1,2,3  (2.26) 

The  demagnetization  factors  Fj  are  given  by 

Fj  =  (/tM  -  1)  /  (1  +  A,  -  firZ]  /  {2^lr2\)  (2.27) 
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Ai  =  A2  -■=  e  (e  +  fJ)  [e}  -  l)“‘  ;  A3  =  -2e  (e"*  -h  E)  (e^  -  l)‘‘  (2.28) 

=  In  -  [e^  -  -  l)  ^  for  e  >  1  (prolate)  (2.29) 

E  —  ^arctan  |l  ~  e^|  ^  for  e  <  1  (oblate)  (2.30) 

Ai  =  A2  =  A3  =  ^  for  c  =  1  (sphere)  (2.31) 


Field  quantities  have  been  derived  in  vhc  body-fixed  coordinate  system  whereas  quanti¬ 
ties  are  needed  in  the  space-fixed  system  in  which  measurements  are  made.  A  body-fixed 
vector  u  and  its  related  space-fixed  vector  U  arc  connected  by  the  relation 

u  =  Xij  (2.32) 

where  A  is  the  well  known  Euler  rotation  tensor  corresponding  to  Euler  angles  (0, 0,Q). 

Equations  2.17  through  2.32  allow  one  to  calculate  the  magnetic  field  B  at  a  point  in 
space  due  to  the  presence  of  the  spheroid,  given  the  size  and  shape  of  the  spheroid  {a,  e), 
the  magnetic  material  propenies  of  the  spheroid  and  the  surrounding  medium  {(iruHr2), 
the  location  of  the  geometric  center  of  the  spheroid  Xo  and  the  orientation  of  the  spheroid’s 
symmetry  axis  with  respect  to  the  space  fixed  system  {$,  4).  However,  the  equations  are 
insensitive  to  (/^ri)/'r2)  provided  >~  100.  Tliis  is  generally  tlie  case  for  ferrous 

materials  in  the  earth’s  field.  For  all  model  calculations  done  in  this  study,  we  have  chosen 
fir\  =  1000,  fir2  —  L  which  are  typical  values  encountered  in  practice. 

In  addition.  Equations  2.20  and  2.26  through  2.32  allow  the  dipole  moment  of  the 
spheioid  in  the  space-fixed  cartesian  frame  M  -  {Mi,  M2,  Ma)^  to  be  calculated  as : 

w  /  {Fi  -t-  [i'3  -  Fi]sin^<?cos^<^}Boi  -f  {[F3  -  ] cost? sin ^ cos (/»}Bo3  \ 

M  =  —  I  {[/<’3  -  /'’i]sin^<?cos<jisin^}Boj  +  {[Fi  -  Bj] cost? sin <? sin (/i}Bo3 

\  (IF,  -  /' i] cost? Sint? cos -f  {F,  +  [F3  -  F,] / 

(2.33) 
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3.  Experimental  Method 


3.1  The  Smart  Magnetometer 

To  obtain  magnetic  data  as  a  function  of  position  in  a  plane,  a  novel  instrument  has  been 
developed,  which  has  been  dubbed  the  “smart  magneiometer”.  It  collects  simultaneous 
magnetic  and  position  data  and  estimates  location  and  identity  of  compact,  axially  symmet¬ 
ric,  ferrous  objects.  The  instrument  is  self-contained  and  can  be  used  in  a  person-portable 
or  vehicle-mounted  role. 

The  second  generation  of  the  instrument  was  used  in  this  study.  It  is  shown  in  Fig.  3.1. 

The  major  components  are  a  total  field  magnetometer,  a  tricycle  can  which  holds 
the  magnetometer  sensor  head  and  has  a  front  wheel-mounted  shaft  encoder  serving  as  a 
linear  position  sensor,  a  Motorola  68332  EVB  microcontroller,  an  interface  between  the 
magnetometer  and  position  encoder  signals  and  the  microcontroller  and  a  lap-top  computer. 
Power  to  the  entire  system  can  be  supplied  by  the  magnetometer  battery  pack,  but  for 
convenience  during  the  many  experiments  described  here,  a  DC  power  supply  was  used. 
Interconnections  between  the  electronic  components  is  shown  in  Figure  3.2,  while  details 
of  the  smart  magnetometer  interface  are  shown  in  Figure  3.3. 

A  model  VlOl  self-oscillating  cesium  vapour  magnetometer,  made  by  Scintnex,  Con¬ 
cord,  ON,  Canada,  is  used  to  obtain  the  magnetic  data.  It  has  a  precision  due  to  quantization 
error  for  a  given  sensor  head  orientation  of  ±0.1  nT  over  a  nominal  range  of  20000  nT  to 
100000  nT.  The  accuracy,  limited  by  heading  error,  is  a  function  of  the  relative  orientation 
of  sensor  optical  axis  and  ambient  field  direction.  It  is  no  worse  than  ±0.5  nT. 

The  magnetometer  takes  45  ms  to  make  a  field  measurement  and  provides  a  digital 
BCD  output  that  is  updated  every  89ms.  The  magnetometer  output  is  buffered  and  fed 
to  paiallel  digital  input  ports  of  the  microcontroller.  As  well,  the  two  phase  outputs  from 
the  position  encoder  are  fed  to  the  input  of  an  advanced  timer/countcr  unit  (TPU)  on  the 
microcontroller.  The  pr  sition  encoder  interrupts  the  microcontroller  and  increments  or 
decrements  a  register  2(XX)  times  for  each  revolution  of  the  front  cart  wheel.  (The  TPU 
compares  the  relative  phases  of  the  encoder  signals  to  determine  the  direction  of  motion 
of  the  wheel  and  correspondingly  increments  or  decrements  the  register.)  The  wheel 
circumference  of  0.618  m  then  corresponds  to  a  limiting  positional  resolution  of 0.309  mm. 
However,  the  effective  resolution  is  governed  by  the  magnetometer  field  measurement  time, 
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Figure  3.1 

Second  generation  “smart  magnetometer”.  The  tricycle  cart  is  on  the 
left.  The  position  encoder  is  attached  to  the  front  wheel  of  the  can.  The 
magnetometer  consists  of  a  display  box  on  the  right,  a  cylindrical  sensor 
head  which  is  attached  to  the  main  rod  of  the  cart  and  an  amplifier  and 
signal  conditioning  unit  and  battery  pack  which  are  not  visible.  Ihe  lap 
top  computer  is  in  the  center  of  the  photograph,  but  the  microcontroller 
is  hidden  from  view. 
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MC68332  EVB  Magnetometer 
Control  Box 


Smart  Magnetometer  Block  Diagram 

Figure  3.2 

Electrical  connections  between  components  of  the  smart  magnetometer. 
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Smart  Magnetometer  interface 

Figure  3.3 

Smart  magnetometer  interface.  Interface  buffers  signals  between  the 
magnetometer  and  position  encoder  on  the  left  and  the  microcontroller 
on  the  right.  Magnetometer  and  microcontroller  signal  names  are  those 
used  by  their  respective  manufacturers  [13],  [  14]. 
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since  the  magnetic  field  value  is  uncertain  within  this  time  period.  The  effective  positional 
resolution  in  the  direction  of  motion  is  thus  dependent  on  cart  speed  and  is  1.12S  cm  at 
0.25  m/sec  and  1.350  cm  at  0.30  m/sec.  The  software  running  on  the  microcontroller  reads 
the  magnetic  field  ports,  converts  the  values  to  a  magnetic  field  and  then  reads  the  position 
register  whenever  the  magnetometer  outputs  a  new  value.  TTie  magnetic  field  conversion 
time  is  less  than  the  spacing  between  encoder  interrupts  (~  1  msec  at  0.3  m/sec)  and  so  the 
position  value  corresponds  to  the  position  at  which  the  magnetometer  signalled  that  it  had 
acquired  the  magnetic  value.  The  magnetic  field  value  and  the  current  position  value  are 
stored  if  a  predetermined  distance  has  been  spanned  since  the  last  update. 

The  microcontroller  is  the  heart  of  the  instrument.  It  controls  collection  of  magnetic 
and  poitition  data  and  algorithm  execution.  The  menu-driven  software  has  the  following 
modes: 

•  Free  running  -  Gives  bai'  graph  indication  of  field  stfength.  Allows  the  operator  to 
find  regions  which  have  magnetic  field  values  of  sufficient  strength  and  spatial  extent 
to  warrant  further  investigation. 

•  Locate  -  Instructs  the  operator  when  and  where  to  make  magnetic  field  and  position 
measurements.  Estimates  location  and  moment  of  the  dipole  associated  with  a 
detected  object  from  the  field  measurements  using  the  noniterative  algorithm  of 
Section  2.2. 

•  Identify  -  Identifies  detected  object  from  its  dipole  moment,  estimated  from  the 
“Locate”  mode,  using  the  pattern  classifier  of  Section  2.3. 

•  Design  set  update  -  If  the  identification  was  carried  out  on  a  known  object,  this  mode 
uses  the  results  to  mexlify  or  add  a  feature  vector  to  the  design  set. 

•  Magnetic  map  transfer  -  Allows  magnetic  maps  to  be  transferred  to  the  lap  top 
computer  hard  disk. 

•  Design  set  transfer  -  Allows  design  sets  to  be  transferred  to/from  the  lap  top  computer 
hard  disk. 

The  lap  top  computer  is  initially  used  to  down  load  code  to  the  n.icrocontroller  and  acts 
as  a  file  server  for  magnetic/position  data  and  design  sets  for  experimental  purposes.  In 
an  fielded  instrument,  code  can  reside  permanently  on  the  microcontroller  EPROM  and  the 
lap  top  can  be  replaced  by  a  hand-held  terminal  or  notebook  computer. 

The  location  algorithm  execution  time  is  a  function  of  how  finely  the  field  is  sampled 
in  the  plane  above  the  object  of  interest  because  this  determines  how  many  magnetic  data 
values  must  be  searched  to  find  the  extrema  ptositions.  The  location  algorithm  executed  in 
five  seconds  for  a  170  cm  (~2  cm  .sample  increments)  by  210  cm  (5  cm  sample  increments) 
grid  38(X)  total  data  points  in  the  magnetic  map).  Execution  lime  for  the  identification 
algorithm,  which  is  independent  of  magnetic  map  size,  took  about  30  seconds  for  a  6  object 
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design  set.  These  times  are  still  small  compared  to  the  data  collection  time.  It  took  6 
minutes  to  collect  the  data  for  a  170  cm  (~2  cm  sample  increments)  by  210  cm  (10  cm 
sample  increments)  grid  and  13  minutes  for  a  170  cm  (~2  cm  sample  increments)  by  210 
cm  (5  cm  sample  increments)  grid. 


3.2  Experimental  Layout 

Measurements  were  made  in  the  Threat  Detection  Group’s  nonmetalHc  laboratory  which 
has  been  described  in  detail  in  [7J.  The  laboratory  is  a  12  m  diameter  hemispherical 
nonmetalHc  building  with  minimal  metal  content.  Magnetic  Held  gradients  are  typically  no 
greater  than  1  or  2  nT/m  in  the  usable  portions  of  the  building.  Magnetic  Held  fluctuations 
are  limited  by  geomagnetic  noise  and  were  typically  ±1-2  nT/hour  during  the  time  of  day 
in  which  measurements  were  made  (mid- morning,  early  afternoon).  In  the  time  necessary 
to  measure  a  magnetic  field  map  (~15  minutes  or  less),  background  noise  fluctuations  were 
less  than  the  quantization  error  of  magnetic  data  received  by  the  microcontroller  (1  nT)  and 
thus  no  reference  magnetometer  was  used  for  background  subtraction.  The  experimental 
set-up  is  shown  schematically  in  Fig  3.4. 

The  measurement  table  on  which  the  magnetometer  moves  was  situated  in  a  low  gradient 
part  of  the  laboratory.  The  table  was  made  of  wood  with  brass  screws  and  was  heavily 
reinforced  to  ensure  negligible  displacement  of  the  measurement  surface  when  walked  on. 
The  horizontal  measurement  surface  was  3.66  m  on  a  side  and  was  1 .04  m  above  the  floor. 
The  table  surface  was  painted  in  alternating  light  and  dark  coloured  guide  lines  spaced  4.8 
cm  apart.  The  table  was  oriented  so  that  the  guide  lines  ran  along  a  direction  5®  from  a 
north-south  direction  (Fig.  3.4).  This  direction  was  labelled  X\.  The  horizontal  direction 
orthogonal  to  the  X\  direction  was  labelled  Xt.  A  start  line,  parallel  to  the  Xi  direction, 
was  situated  at  either  end  of  the  table.  Because  of  practical  considerations  such  as  the  finite 
size  of  the  magnetometer  cart,  placement  of  the  start  lines  and  positioning  of  the  operator, 
the  usable  region  of  the  measurement  surface  was  about  2.7m  in  the  X'l  direction  and  2.55m 
in  the  X2  direction.  The  intersection  of  the  lower  start  line  and  the  left-most  guide  line 
was  considered  to  be  (0,0,  -0.24m)  in  the  space-fixed  coordinate  system  since  the  nominal 
center  of  the  magnetometer  active  volume  was  24  cm  above  the  table  surface. 

The  object  of  interest  was  placed  in  a  wooden  holder  that  allowed  the  object  to  be  rotated 
independently  about  a  vertical  or  horizontal  axis  and  locked  at  a  desired  orientation.  The 
geometric  center  of  the  object  of  interest  was  then  placed  directly  below  the  center  of  the 
measurement  table  at  a  variable  depth  and  orientation.  The  orientation  of  the  object  was 
defined  by  two  angles.  ITie  polar  angle,  0,  was  the  angle  between  the  X3  (vertical)  direction 
and  the  symmetry  axis  of  the  object.  The  azimuthal  angle,  was  the  angle  between  the 
north-south  direction  and  the  projection  of  the  symmetry  axis  on  the  X1X2  (horizontal) 
plane. 

A  separate  wooden  table,  on  which  the  electronics  was  placed,  was  situated  roughly  2 
in  from  the  measurement  table.  Although  all  of  the  electronics  could  have  been  carried  by 
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start  tins 
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display  box 


ELECTRONICS 

TABLE 

power  supply' 


Interface  card 

Figure  3.4 


microcontroller 


Top  view  of  the  experimental  set-up  for  measurement  of  magnetic  fields 
of  compact  objects. 
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the  operator,  it  was  felt  better  for  experimental  purposes  to  keep  the  components  separate 
for  ease  of  modification  and  to  keep  them  distant  from  the  measurement  surface  to  ensure 
that  stray  magnetic  fields  from  the  electronics  would  not  interfere  with  the  magnetometer 
sensor.  Because  of  cable  length  restrictions,  the  magnetometer  electronics  box  was  kept  on 
the  measurement  table  approximately  0.5  m  from  the  magnetometer  cart. 


3.3  Procedure 


To  collect  a  magnetic  map  the  following  procedure  was  used.  Prior  to  collecting  the 
initial  map,  the  electronics  were  connected,  powered  up  and  allowed  to  stabilize  for  at  least 
ten  minutes.  The  operator  removed  all  metal  and  placed  it  far  from  the  measurement  table. 
The  microcontroller  code  was  downloaded  from  the  lap  top  computer  and  the  program  was 
initiated.  For  all  maps,  the  object  of  interest  was  placed  in  the  appropriate  position  and 
orientation.  The  magnetometer  cart  was  positioned  with  its  leai’  wheels  at  the  lower  start 
line  and  its  front  wheel  on  the  lefLmost  guide  line  of  Figure  3.4  (cart  pointing  in  the  positive 
X\  direction).  The  “Locate”  mode  of  the  software  was  selected  and  initial  parameters  were 
entered  (vVi,  Xj  maximum  dimensions,  sample  spacing).  The  background  magnetic  field 
was  sampled  for  approximately  2  seconds  and  then  a  tone  signalled  the  operator  to  roll  the 
cart  with  the  front  wheel  running  along  the  guide  line.  Data  would  not  be  collected  until 
the  wheels  started  to  roll.  While  rolling,  simultaneous  magnetic  and  position  information 
were  automatically  acquired.  When  the  microcontroller  sensed  that  sufficient  distance  had 
been  traversed,  data  collection  was  halted  and  the  operator  was  signalled  to  stop.  The  set 
of  position  and  magnetic  data  collected  along  a  guide  line  is  called  a  “scan”  or  a  “scan 
line”.  The  operator  then  set  the  cart  with  its  rear  wheels  at  the  upper  start  line  and  its  front 
wheel  on  the  next  guide  line  in  the  positive  A'2  direction  (cart  pointing  in  the  negative 
direction).  The  software  had  a  built-in  delay  during  which  data  collection  was  suspended 
to  allow  such  positioning  of  the  cart.  The  microcontroller  would  then  signal  the  operator 
to  commence  rolling  the  cart.  (Note  that  the  microcontroller  software  can  use  the  wheel 
position  encoder  to  guide  the  magnetometer  to  the  start  of  the  next  scan.  This  is  not 
as  accurate  as  using  well  defined  guidelines  and  starting  points  and  so  this  feature  was  not 
used  for  these  experiments.)  The  procedure  was  repeated  and  a  two  dimensional  anay  of 

iiiE^ilCuC  VaiUcS  vciSuS  y\2  vviti  vK'  vaitvU  o  map/  Wao  ill  a 

raster  scan  fashion  until  the  microcontroller  signalled  the  operator  to  halt.  The  location 
of  the  object  and  its  associated  dipole  moment  were  then  estimated  and  the  object  was 
classified.  If  desired,  magnetic  data  were  transfeiTed  to  hard  disk  for  permanent  storage 
and  then  the  prt.cedure  was  repeated  for  a  new  map. 


3.4  Initial  Calibration 

Prior  to  commencing  the  location  and  identification  experiments,  a  number  of  magnetic 
maps  were  collected  to  optimize  operational  procedures  and  parameters.  Yawing  of  the 
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se  nsor  head  was  initially  noted  as  a  problem,  as  was  an  offset  between  the  magnetic  sensor 
head  path  and  the  encoder  wheel  path.  Both  of  these  were  remedied  by  modifying  the  cart. 

Optimum  scanning  speed  for  the  sensor  was  determined.  Because  of  the  fixed  temporal 
spacing  between  magnetometer  outputs,  the  cart  speed  had  to  be  adjusted  so  that  the  sample 
spacing  was  not  too  large,  while  still  quickly  completing  a  scan.  A  sample  spacing  of  about 
2  cm  was  desired.  It  was  found  that  a  reasonable  sensor  speed  was  roughly  0.30  to  0.25 
m/sec.  llie  cart  would  require  8  to  10  seconds  to  move  from  start  to  finish  of  one  scan  line 
and  the  microcontroller  would  store  a  magnetic  field  reading  every  2.2  to  2.7  cm  along  the 
Xi  direction. 

A  magnetic  map  was  obtained  with  no  object  in  place,  in  order  to  measure  the  back¬ 
ground  magnetic  field.  Tlie  map  is  shown  in  Figure  3.5.  Hie  field  gradient  is  seen  to 
be  roughly  constant  and  can  be  characterized  by  a  single  component  pointing  from  the 
northeast  to  the  southwest  corner  of  the  table  with  magnitude  ~  2  nT/m.  The  effect  of 
the  background  is  small  but  significant,  particularly  for  the  moment  estimation  when  the 
magnetic  field  of  the  object  is  weak,  lypical  examples  of  the  difference  in  estimation  with 
and  without  background  subtraction  are  shown  in  Tables  I  and  II.  The  object  was  spheroid 
“F”  (see  Table  IV  for  description)  in  a  vertical  orientation  at  depths  of  0.598  m  (Table  I) 
and  0.802  m  (Table  \I).  Peak  field  value  due  to  the  object  was  ~  lOOOnT  for  Table  I  and 
~  350nT  for  Table  II.  For  both  Tables,  the  sample  spacing  in  X2  direction  was  4.8  cm. 
Seven  field  maps  were  used  for  estimation.  Standard  deviation  for  the  magnitude  is  actually 
the  root  average  square  distance  (intraset  distance)  between  vectors,  based  on  7  field  maps 
for  Table  I  and  3  field  maps  for  Table  II.  The  difference  in  estimates  between  corrected 
and  uncorrected  maps  is  at  most  only  slightly  greater  than  can  be  accountwl  for  by  the 
standard  deviations  for  the  0.598  m  depth.  The  difference  becomes  significantly  greater 
than  the  statistical  uncertainty  at  a  depth  of  0.802  m,  which  is  to  be  expected  because  of  the 
diminished  field  due  to  the  object. 

In  the  experiments  that  follow,  the  background  map  was  subtracted  from  the  magnetic 
field  map  due  to  the  object  before  applying  the  location  and  identification  algorithms. 

The  measurement  table  guide  lines  would  allow  scan  lines  a ;  close  together  as  4.8  cm 
or  at  any  multiple  of  4.8  cm.  It  was  desirable  to  use  the  closest  spacing  possible  for  the 
study  of  the  location  algorithm,  but  it  was  necessary  to  know  how  sensitive  the  algorithm 
estimates  were  to  the  X2  spacing.  To  test  this  sensitivity,  a  set  of  measurements  were  made 
on  a  fixed  object  at  a  fixed  oiientation.  The  depth  chosen  was  small  (0.'^9  m)  since  the 
estimation  should  be  more  sensitive  to  grid  spacing  when  field  maps  aie  narrow  as  they  are 
at  shallow  depths.  The  results  are  shown  in  Tablelll,  The  object  chosen  was  spheroid  “F” 
(see  Table  IV)  alO  --  90°,  </>  =  175°  and  at  a  depth  of  0.39  m.  Five  field  maps  were  used 
for  estimation  of  5  cm  spacing  data  and  two  field  maps  were  u.sed  for  estimation  of  the  10 
cm  spacing  data.  (These  data  were  obtained  on  an  older,  smaller  measurement  table  which 
had  a  minimum  scan  line  spacing  of  5  cm,  rather  than  4.8  cm.)  The  standard  deviation  for 
the  magnitude  is  actually  the  root  average  square  distance.  Clearly,  there  is  not  a  significant 
difference  between  the  two  spacings. 
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Background  on  Big  Table  16/07/92  PM 


Figure  3.5 

Background  magnetic  map  collected  by  the  smart  magnetometer  on  the 
measurement  table  surface. 
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Component 

Uncorrected  Map 
Mean  |  Std.  Deviation 

Corrected  Map 

Mean  |  Std.  Deviation 

Location  Vector  Estimate  (m) 

1 

1.358 

0.010 

1.350 

0.010 

2 

1.277 

0.001 

1.278 

0.002 

3 

-0.540 

0.001 

-0.554 

0.002 

magnitude 

1.941 

0.010 

1.940 

0.011 

Moment  Vector  Estimate  (A  m^) 

1 

0.139 

0.035 

0.185 

0.041 

2 

0-005 

0.003 

0.003 

0.007 

3 

0.845 

0.008 

0.909 

0.0188 

magnitude 

0.857 

0.036 

0.927 

0.046 

Table  1 

Comparison  of  location  and  identification  estimation  results  with  and 
without  background  correction.  Object  was  spheroid  “F”  in  a  vertical 
orientation  at  a  depth  of  0.598  m. 


Component 

Uncorrected  Map 
Mean  [  Std.  Deviation 

Corrected  Map 

Mean  |  Std.  Deviation 

1  Location  Vector  Estimate  (m) 

1 

1.383 

0.006 

1.362 

0.006 

2 

1.322 

0.001 

1.333 

0.001 

3 

-0.708 

0.(X)4 

-0.746 

0.007 

magnitude 

2.040 

0.007 

2.047 

0.009 

Moment  Vector  Estimate  (A  m^) 

1 

0  077 

noni 

A 

0.173 

0.023 

2 

0.008 

0.003 

0.013 

0.005 

3 

0.774 

0.012 

0.893 

0.022 

magnitude 

0.778 

0,012 

0.910 

0.033 

Table  II 

Comparison  of  location  and  identification  estimation  results  with  and 
without  background  correction.  Object  was  spheroid  “F”  in  a  vertical 
oricniation  at  a  depth  of  0.802  m. 
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Component 

Sem  Map 

Mean  |  Std.  Eleviadon 

10  cm  Map 

Mean  |  Std.  Deviation 

1  Location  Vector  Estimate  (m) 

1 

0.865 

0.007 

0.865 

0.012 

2 

1.058 

0.009 

1.066 

0.019 

3 

-0.415 

0.007 

-0.422 

0.027 

magnitude 

1.428 

0.019 

1.437 

0.035 

Moment  Vector  Estimate  (A  m^) 

1 

0.266 

0.020 

0.275 

0.033 

2 

-0.004 

0.017 

-0.002 

0.038 

3 

0.215 

0.011 

0.216 

0.028 

magnitude 

0.342 

0.028 

0.349 

0.058 

Table  III 

Comparison  of  location  and  identification  estimation  results  for  magnetic 
maps  collected  with  an  Xi  spacing  between  scans  of  5  cm  and  10  cm. 
Object  was  spheroid  “F’  at  <?  =  90”  and  (j>  —  175®  at  a  depth  of  0.39  m. 
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4.  Experimental  Results 


Using  the  experimental  method  of  the  last  Chapter,  magnetic  maps  were  obtained  for 
two  mild  steel  solid  spheroids  for  a  variety  of  depths  and  orientations.  The  two  objects 
are  summarized  in  Table  IV.  Both  objects  arc  typical  of  the  size  and  shape  of  unexploded 
ordnance.  In  the  following  discussion,  the  type  of  spheroid  and  its  orientation  for  a  paiticular 
experiment  will  be  designated  by  a  letter  (assigned  for  historical  reasons)  and  two  integers. 
For  example,  M  90  175  refers  to  the  type  “M”  spheroid  with  a  polar  angle  of  90“  and 
azimuthal  angle  175“.  The  orientation  for  a  vertical  spheroid  is  deni.ied  0  0,  even  though 
the  azimuthal  angle  is,  strictly  speaking,  undefined.  The  orientation  angles  are  defined  in 
Section  3.2. 

Figures  4.1  to  4.4  are  typical  of  the  variation  of  magnetic  field  maps  for  a  single  object 
for  different  orientations.  Both  the  maximum  field  value  and  the  relative  position  of  the 
maximum  and  minimum  field  value  change  with  orientation.  The  maximum  field  value  is 
seen  to  decrease  with  increasing  polar  angle.  This  illustrates  that  simple  location  schemes 
which  assume  that  the  object  center  is  directly  under  the  maximum  field  position  will  have 
serious  errors  as  will  simple  identification  schemes  which  use  only  the  maximum  field  value 
as  a  feature. 


The  location  estimates  for  all  experiments  are  given  in  Table  V.  The  dipole  moment 
estimates  for  all  experintents  are  given  in  Table  VI.  For  both  tables,  where  multiple 
experiments  were  done  for  identical  conditions,  the  standard  deviations  are  given.  Note 
that  the  range  of  depths  is  from  0.39  m  to  0.802  m.  The  lower  limit  was  mainly  impxrscd  by 
the  measurement  table  surface  thickness  plus  magnetometer  height  plus  minimum  object 
dimension  (  ~  4  +  24  +  10,5  cm  respectively),  while  the  uppjer  limit  was  set  by  the 
requiremeiit  to  havea  sufficient  signal-io-noise  ratio  (SNR)  to  get  reasonable  precision  for 
the  spheroids  studied. 


It  is  necessary  to  characterize  both  the  precision  and  the  accuracy  of  the  location 
algorithm,  lb  determine  the  precision  or  repeatability,  one  can  look  at  the  location  and 
diprolc  moment  estimates  for  independent  magnetic  maps  that  were  obtained  for  the  same 
object,  orientation  and  depth.  Precision  of  a  location  vector  compxrnent  was  measured  by  the 
standard  deviation  of  the  comp)onent,  expressed  as  a  fraction  of  the  depth,  while  precision 
of  a  moment  vector  compxjiient  was  measured  by  the  standard  deviation  of  the  compronent, 
expres.sed  as  a  percentage  of  the  average  magnitude  of  the  vector.  A  good  measure  of  the 
deviation  of  vectors  sampled  from  the  same  population  is  the  intraset  distance.  Given  a  .set 
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Object  Label 

Type 

a(m) 

e 

F 

Prolate  Spheroid 

0.0525 

3.5 

M 

Prolate  Spheroid 

0.05896 

2.0 

Table  IV 

Two  objects  used  for  experiments.  Spheroid  dimensions  are  defined  in 
Figure  2.2 


Prolote  Spheroid  TOO  d=--0.59arn  14/05/92  PM 


Measured  total  magnetostatic  field  map  versus  position  in  a  horizontal 
plane  for  spheroid  F  at  <?  =  0°.  Depth  is  0.598  ni.  Ambient  (earth’s)  field 
magnitude  has  been  subtracted. 
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Figure  4.2 

Measured  total  magnetostatic  field  map  vcr-”s _ _  .*  u  • 

plane  for  spheroid  F  at  0  ~  30"  anH  a  *  "O^^ontaJ 

Ara6fcn.(eanh’.)  field  magninicte  ha.  be.n'subttacted'!^'’''’  " 
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Prolate  Spheroid  F  60  175  d=0.525m  06/08/92  PM 


Measured  total  magnetostatic  Held  map  versus  position  in  a  horizontal 
plane  for  spheroid  F  at  5  =  60“  and  (f>  ==  175®.  Depth  is  0.525  m. 
Ambient  (earth’s)  field  magnitude  has  been  subtracted. 
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Prolote  Spheroid  F  90  175  d-0.655m  16/07/92  PM 

jul  6_2.dat 


Measured  total  magnetostatic  Held  map  versus  position  in  a  hori  cental 
plane  for  spheroid  F  at  <1  =  90°  and  (j>  =  175°.  Depth  is  0.655  m. 
Ambient  (earth’s)  field  magnitude  has  been  subtracted. 
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spheroid 

TVpe. 

Orientaiion 


M  10  175 


F  90  175 


F  60  175 


F  30  175 


F  30  265 


M  90  265 


Measured 
■^01  A'o2  A'o3 

(m)  (m)  (m) 


0.87 


1.321 


1.357 


1.357 


1.347 


0.87 


1.362 


0.855 


0.87 


1.36 


1.362 


1.362 


1.362 


1.362 


0.87 


0.87 


1.04 


1.288 


1.333 


1.333 


1.293 


1.04 


1.313 


1.04 


1.04 


1.33 


1.313 


1.313 


1.313 


1.313 


1.04 


1.04 


0.5G 


0.598 


0.698 


0.802 


0.565 


0.50 


A'oi 

(m) 

Estimated 

a:o2 

(m) 

ATos 

(m) 

0.870 

1.05 

0.411 

1.350  ±0.010 

1.278  ±  0.002 

0.554  ±  0.002 

1.372  ±0.002 

mwmamrM 

0.642  ±  0.000 

1.362  ±0.006 

1.3.33  ±0.001 

0.746  ±  0.006 

1.344  ±0.002 

1.298  ±0.000 

0.487  ±  0.001 

0.862 

1.05 

0.458 

1.368 

1.311 

0.534 

0.865  ±0.007 


0.894 


1.391 


1.392 


1.050 


1.306 


1.307 


1.301 


1.389 


0.884  ±0.005 


0.877 


1.326 


1.040  ±0.001 


1.070 


1.003 


0.504 


0.676 


0.539 


0.494 


0.556 


0.499 


0.453  ±  0.007 


0.455 


0.524 


Table  V 

Location  estimation  results  using  experimental  data  from  spheroids. 
Where  multiple  experiments  have  been  done  under  identical  conditions, 
standard  deviations  are  given. 
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Spheroid 

Type, 

Orientation 

depth 

(m) 

Mathc 

Ml 

(Am^) 

:matical  ^ 
M2 

(Am^) 

dodel 

M3 

(A-m^) 

Ml 

(A-m^) 

Estimated 

M2 

(A-m^) 

Ms 

(Am^) 

FOO 

0.500 

0.065 

0.000 

1.069 

0.072 

0.001 

0.716 

0.698 

0.065 

0.000 

1.069 

0.177 

±0.005 

0.013 
±  0.013 

0.876 

±0.001 

0.565 

0.065 

0.000 

1.069 

0.161 
±  0.003 

-0.009 

±0.001 

0.701 

±0.002 

0.598 

0.065 

0.000 

1.069 

0.185 

±0.042 

0.003 

±0.007 

0.909 
±  0.019 

0.802 

0.065 

0.000 

1.069 

0.173 
±  0.023 

0.013 

±0.005 

0.893 
±  0.022 

F  90  175 

0.390 

0.325 

-0.023 

0.212 

0.266 
±  0.020 

-0.004 
±  0.017 

0.215 
±  0.01 1 

0.440 

0.325 

-0.023 

0.212 

0.315 

-0.001 

0.245 

0.655 

0.325 

-0.023"^ 

0.212 

0.238 

-0.015 

0.182 

0.525 

0.325 

0.023 

0.212 

0.410 

-0.035 

0.191 

M  90  265 

0.438 

0.059 

0.007 

0.189 

0.082 

±0.006 

-0.061 

±0.001 

0.158 

±0.007 

0.468 

0.059 

0.007 

0.189 

0.039 

0.029 

0.105 

0.580 

0.059 

0.007 

0.189 

0.028 

-OO^ 

0.123 

F  60  175 

0.530 

0.621 

-0.049 

0.539 

0.779 

-0.062 

0.550 

F  30  175 

0.525 

0.499 

-0.038 

0.967 

0.554 

-0.040 

0.833 

F  60  265 

0.525 

0.099 

0.386 

0.436 

0.137 

0.495 

0.541 

F  30  265 

0.525 

0.098 

0.375 

0.864 

0.169 

0.028 

0.817 

M  10  175 

0.500 

0.105 

-0.004 

0.455 

0.030 

-0.003 

0.300 

Tabic  VI 

Dipole  moment  estimation  results  using  experimental  data  from 
spheroids.  Where  multiple  experiments  have  been  done  under  identical 
conditions,  standard  deviations  are  given.  Theoretical  moment  compo¬ 
nent  values  are  obtained  from  the  mathematical  model  of  Section  2.4. 
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of  vectors,  a\i  —  1, A',  the  intraset  distance,  'D^,  is  the  mean  square  distance  between 
vectors  given  by  [15] 


= 


1 

A'(A-l) 


i=i  k=i 


(4.1) 


where  a\  is  the  A;th  component  of  a*  and  k  —  The  precision  in  estimation  of  location 

vectors  was  taken  to  be  the  square  root  of  the  intraset  distance  between  individual  location 
vector  estimates  for  a  given  object,  orientation  and  depth.  For  two  depths,  O.SOO  m  and 
0.S2S  m,  different  orientations  and/or  objects  were  used  to  determine  the  precision  in  depth 
whereas  for  other  depths,  these  parameters  were  constant  for  a  given  depth.  The  precision 
in  estimation  of  moment  vectors  was  assigned  to  be  the  square  root  of  the  intraset  distance 
between  individual  moment  vector  estimates.  Since  moments  are  a  function  of  object  type 
and  orientation,  precision  in  moments  could  not  be  determined  for  the  O.SOO  and  0.525  m 
depths.  The  results  of  the  precision  in  estimation  are  given  in  Tables  VII  and  VIII  and 
Figures  4.5  and  4.6. 


From  Tables  Vlland  VIII  and  Figures  4.5  and  4.6,  it  can  be  seen  that  the  precision  in 
estimating  location  varies  from  0.12  to  1.90%.  However,  there  is  also  a  variation  in  location 
estimate  with  object  type  and/or  orientation,  since  the  two  depths  with  the  worst  precision 
(1.67,  1.90%)  are  the  only  ones  which  have  different  types  and  orientations  included  in 
the  precision  estimates.  If  these  are  excluded,  the  precision  in  depth  ranges  from  0.12  to 
1.35%.  The  precision  in  estimating  the  dipole  moment  varies  from  0.51  to  8.21  %,  There 
appears  to  be  a  slight  increase  in  uncertainty  as  depth  decreases,  but  the  precision  is  not 
closely  correlated  with  the  peak  magnetic  held,  litis  variation  in  precision  may  be  partly 
because  the  widths  of  the  magnetic  field  maps  scale  linearly  with  depth.  When  sampling 
with  a  fixed  width  sampling  grid,  estimates  will  be  more  sensitive  to  narrow  maps.  Also, 
the  lowest  depth  cases  were  done  for  spheroids  with  a  horizonial  symmetry  axis.  The  shape 
of  the  magnetic  map  for  these  cases  will  differ  from  those  for  tne  remainder  which  arc  for 
vertical  axis  spheroids  (compare  Figs.  4.1  and  4.4).  The  use  of  a  fixed  width  sampling  grid 
may  then  again  make  estimates  more  sensitive  to  one  orientation  than  another. 


It  is  also  interesting  to  note  that  the  precision  of  the  dipole  moment  estimates  is  always 
larger  than  the  precision  of  the  location  estimates  for  the  same  object,  orientation  and  depth. 
Tnis  is  as  expected,  since  Equation  2.S  siiows  that  for  a  given  orientation,  the  magnitude  of 
the  moment  is  equal  to  the  cube  of  the  depth  times  a  constant.  Simple  enor  propagation 
theory  then  suggests  that  the  percentage  error  in  moment  should  be  three  times  that  in  depth. 
Since  error  in  depth  usually  dominates  the  error  in  the  other  directions  [10],  the  percentage 
error  in  dipole  moment  should  be  at  least  three  limes  the  percentage  error  in  location. 


Although  the  precision  (repeatability)  is  quite  good,  it  is  also  necessary  to  know  what  is 
the  accuracy  in  estimation;  that  is,  are  the  estimates  biased?  To  get  a  preliminary  assessment 
of  accuracy,  the  depth  estimates  averaged  over  all  object  types  and  orientations  for  a  given 
depth  were  plotted  versus  measured  depth.  The  results  are  presented  in  Table  IX  a)id 
Figure  4.7 
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Depth 

(ra) 

Spheroid 

TVpe, 

Orientation 

Number 

of 

Trials 

Peak 

Field 

(nT) 

(%) 

Xoa 

(%) 

Xo3 

(%) 

Xo 

(%) 

0.390 

F90  175 

5 

~600 

1.79 

2.31 

1.79 

1.35 

0.438 

M  90  265 

4 

~300 

1.14 

0.23 

1.60 

0.86 

0.500 

F.M-- 

2 

- 

0.78 

0.12 

4.69 

1.67 

0.525 

F-- 

5 

- 

4.27 

1.95 

5.43 

1.90 

0.565 

FOO 

2 

~1000 

0.35 

0.00 

0.17 

0.12 

0.598 

FOO 

7 

~800 

1.67 

0.33 

0.33 

0.56 

0.698 

FOO 

2 

-600 

0.28 

0.28 

~a6^ 

0.17 

0.802 

FOO 

3 

-350 

0.75 

0.12 

0.75 

0.44 

Table  VII 

F*recision  in  location  vector  estimation  for  the  location  algorithm  of  Sec¬ 
tion  2.2  using  experimental  data.  For  a  given  depth,  the  precision  in  a 
component  Xoi,i  -=  1,2, 3  is  the  standard  deviation  for  the  component 
expressed  as  a  percentage  of  the  depth.  Precision  of  the  vector  Xq  is 
the  square  root  of  the  intraset  distance  for  all  vectors  at  a  given  depth 
expressed  as  a  percentage  of  the  vector  magnitude.  Data  at  0.500  and 
0.525  m  depths  includes  different  orientations  and/or  objects. 
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Depth 

(m) 

Spheroid 

T^'pe, 

Orientation 

Number 

of 

Trials 

Peak 

Field 

(nT) 

Ml 

(%) 

Ml 

(%) 

Ms 

(%) 

M 

(%) 

0.390 

F  90  175 

5 

-600 

5.85 

4.97 

3.22 

8.21 

0.438 

M  90  265 

4 

-300 

3,20 

0.45 

3.73 

6.65 

0.565 

FOO 

2 

-1000 

0.40 

0.13 

0,27 

0.51 

0.598 

FOO 

7 

-800 

4.54 

0.75 

2.04 

4.97 

0.698 

FOO 

2 

-600 

0.56 

1.46 

0.11 

1.58 

0.802 

FOO 

3 

-350 

2.52 

0.53 

2.41 

2.61 

Table  VIII 

Precision  in  dipole  moment  vector  estimation  for  the  location  algorithm 
of  Section  2.2  using  experimental  data.  For  a  given  depth,  the  precision  in 
a  component  M,  ,  i  =  1, 2, 3  is  the  standard  deviation  for  the  component 
expre:sed  as  a  percentage  of  the  average  moment  magnitude.  Precision 
of  the  vector  is  the  the  square  root  of  the  intraset  distance  for  all  vectors 
of  a  given  object  type,  orientation  and  depth,  expressed  as  a  percentage 
of  the  vector  magnitude. 


Measured  Depth 
(m) 

Estimated  Depth 
(in) 

Measured  Depth 
(m) 

Estimated  Depth 
(m) 

0.390 

0.415  ±0.007 

0.565 

0.487  ±0.001  1 

0.438 

0.453  ±  0.007 

0.580 

0.524 

0.440 

0.504 

0.598 

0.554  ±  0.002 

0.468 

0.455 

0.655 

0.676 

r\  trr\f\ 

U.JUU 

0.435  ±  0.023 

0.698 

0.642  ±  O.OCX) 

0.525 

0.525  ±  0.029 

0.802 

0.746  ±  0.007 

Table  IX 

Spheroid  depth  estimates  as  a  function  of  depth  for  experimental  data. 
All  objects  and  orientations  are  grouped  together  for  a  given  depth.  Un¬ 
certainties  (precision)  arc  given  if  available. 
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mcusurcd  depth  (m) 


Figure  4.5 

Precision  in  location  vector  estimation  for  the  location  algorithm  using 
experimental  data.  The  precision  is  the  root  mean  square  error  for  all 
vectors  at  a  given  depth.  Data  at  0.500  and  0.525  m  depths  includes 
different  orientations  and/or  objects. 
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measured  depth  (m) 

Figure  4.6 

Precision  in  dipole  moment  vector  estimation  for  the  location  algorithm 
using  experimental  data,  llie  precision  is  the  root  mci  ’  square  eiror  for 
all  vectors  of  a  given  object  type  and  orientation  and  depi.h  expressed  as 
a  fraction  of  the  dipole  moment  vector  magnitude. 
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mcosured  depth  (m) 


Figure  4.7 

Spheroid  depth  estimates  us  a  function  of  depth  for  cxpci'iuicniiti  data. 
All  objects  and  orientations  arc  grouped  together  for  a  given  depth.  Un¬ 
certainties  are  shown  where  available  and  are  in  most  cases  smaller  than 
the  data  points.  The  dashed  line  is  the  line  along  which  data  points  would 
lie  if  the  estimates  were  unbiased. 


DRES-SR-582 


UNCLASSIFIED 


40 


UNCLASSIFIED 


From  the  data,  it  can  be  concluded  that  there  is  a  depth  dependent  bias  in  the  estimation 
of  depth.  The  bias  is  larger  than  can  be  explained  by  statistical  fluctuations.  In  general, 
depth  estimates  are  slightly  high  for  small  depths  and  slightly  low  for  large  depths.  If  one 
examines  Table  V  more  closely,  it  can  be  seen  that  the  algoiithm  generally  underestimates 
the  depth  for  0  <  30“  and  overestimates  for  0  >  60“ .  An  excepdon  is  the  (90  265) 
orientation  for  which  one  depth  is  underestimated  and  two  are  overestimated.  It  is  also 
interesting  to  note  that  the  average  error  in  depth  is  larger  for  vertical  spheroids  than  it 
is  for  higher  polar  angles  (Table  V).  One  explanation  is  that  the  negative  peak  is  smaller 
in  magnitude  and  broader  for  the  vertical  case  than  for  larger  polar  angles,  which  might 
make  the  error  in  csdmating  the  peak  position  laiger.  This  can  be  discounted  since,  if  it 
were  true,  the  precision  in  estimates  for  the  vertical  spheroid  would  be  worse  than  for  the 
horizontal  spheroid.  However,  the  precision  for  the  two  orientations  are  similar.  A  more 
likely  explanation  is  that  for  a  given  orientation,  the  octupole  field  causes  the  centroids  of 
the  magnetic  field  to  be  shifted  from  the  position  of  the  field  extrema,  lliis  shift  will  be 
a  function  of  the  ratio  of  octupole  to  dipole  field  values  in  the  measurement  plane,  which 
will  in  turn  be  dependent  on  the  object  orientation.  Likewise,  the  effect  of  a  small  amount 
of  remnant  magnetization  will  also  be  dependent  on  orientation  [7]. 

To  examine  the  bias  in  estimation  further,  it  is  insti  uctivc  to  look  at  the  bias  in  estimation 
of  the  individual  components  of  the  location  vector.  TTiesc  are  presented  in  Figures  4.8 
to  4. 10 

It  is  seen  that  all  three  components  have  a  depth  dependent  bias  which  is  larger  than 
can  be  accounted  for  by  the  precision  of  the  algorithm,  llie  Xoi  component  of  the  location 
is  almost  always  overestimated.  The  Xi)2  component  is  generally  overestimated  for  depths 
less  than  0.55  m  and  underestimated  for  depths  between  0.55  m  and  0.8m.  Over  the  limited 
depth  range,  the  data  is  consistent  with  a  bias  which  is  periodic  with  depth,  having  a  0.25 
m  pericxi.  Further  experiments  at  larger  and  smaller  depths  might  confirm  this,  A  trend 
similar  to  that  of  the  X^i  component  appears  to  occur  for  the  A'^03  component  except  that 
the  crossover  from  overestimation  to  underestimation  occurs  at  about  0.45  m  depth.  'Iltc 
data  may  be  consistent  with  a  bias  which  is  peruxiic  with  depth,  but  the  uncertainties  in 
estimation  make  deiennination  of  the  peritxi  difficult.  Errors  in  A03  are  generally  hirger 
than  those  for  A"oi  and  A02.  'lire  root  mean  square  (nns)  error  over  all  depths  is  0.029  m 
for  Aoi.  0.019  m  for  Ao2  and  0.045  m  for  A'o3-  This  should  be  compared  with  the  average 
precision  for  a  fixed  object  and  orientation  which  is  0.005  m  for  A'oi,  0.003  m  for  A02,  and 
0.004  m  for  A'^us- 

Analy.sing  the  bias  in  the  estimation  of  the  dipole  moments  is  more  difficult,  since  the 
“true”  moment  components  cannot  be  directly  measured  as  can  the  location  components. 
We  can,  however,  compare  the  moment  component  estimates  for  a  fixed  object  type  and 
orientation  with  theoretical  values  obtained  from  the  model  of  Section  2.4.  Figures  4.11 
to  4.13  show  the  variation  with  depth  of  the  estimates  of  the  three  moment  comimnents 
for  a  type-F  spheroid  with  orientation  0  0.  Figures  4.14  to  4.16  show  the  same  thing  for  a 
type  F  spheroid  with  orientation  90  175  and  Figures  4.17  to  4.19  show  the  variation  for  a 
type-M  spheroid  with  orientation  90  265.  (Also  shown  on  each  Figure  are  estimates  of  the 
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measured  depth  (m) 


Figure  4.8 

Error  in  estimates  of  A'oi  component  as  a  function  of  spheroid  depth  for 
experimentai  data.  Error  is  denned  as  estimated  A'm  minus  measured 
Aoi .  All  objects  and  orientations  are  grouped  together  for  a  given  depdt. 
Uncertainties  are  shown  where  available  and  are  in  some  cases  smaller 
than  the  data  points.  The  dashed  line  is  the  line  along  which  data  points 
would  lie  if  the  estimates  were  unbiased. 
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Figure  4.9 


Error  in  estimates  of  A'^02  component  as  a  function  of  spheroid  depth 
for  experimental  data.  Error  is  defined  as  estimated  component  minus 
measured  component.  All  objects  and  orientations  are  grouped  together 
f(.r  a  given  depth.  Uncenainiies  are  .shown  where  available  and  are  in 
some  cases  smaller  than  the  data  points.  The  dashed  line  is  the  line  along 
which  data  points  would  lie  if  the  estimates  were  unbiased. 
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Figure  4.10 

Error  in  estimates  of  A''o3  component  as  a  function  of  spheroid  depth 
for  experimental  data.  Error  is  defined  as  estimated  component  minus 
measured  component.  Alt  objects  and  orientations  are  grouped  together 
for  a  given  depth.  Uncertainties  are  shown  where  available  and  are  in 
some  cases  smaller  than  the  data  points.  The  dashed  line  is  the  line  along 
which  data  points  would  lie  if  the  estimates  were  unbiased. 
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dipole  moment  components  using  as  input  magnetic  field  maps  derived  from  the  model  of 
Section  2.4.  Discussion  of  these  will  be  deferred  to  the  next  Chapter.) 

A  number  of  points  should  be  noted  about  Figures  4.1 1  to  4.19.  ITiere  is  a  deviation 
between  theoretical  values  and  experimental  estimates  that  is  greater  than  the  uncertainties  in 
the  estimates.  Within  the  confines  of  the  limited  number  of  data  points,  the  deviation  does  not 
exhibit  a  clear  trend  with  depth.  To  explore  a  possible  trend  with  object  type  or  orientation, 
the  difference  between  estimated  and  theoretical  moment  components,  normalized  by  the 
magnitude  of  the  theoretical  moment,  were  plotted  for  various  objects  and  orientations. 
For  multiple  estimates  of  the  same  object  type  and  orientation  at  different  depths,  the 
unweighted  average  of  estimates  at  all  depths  for  that  object  type  and  orientation  was  used 
and  the  uncertainty  in  estimate  was  taken  to  be  the  standard  deviation  of  the  estimates  about 
the  average.  These  results  are  shown  in  Figures  4.20  to  4.23. 

The  unweighted  average  difference  between  estimated  and  theoretical  components  as  a 
percentage  of  the  theoretical  dipole  moment  is  4.8  ±  7.6%  for  Mj,  —3.3  ±  15.5%  for  M2, 
—  10.8  ±  17.5%  for  Mi  and  24.5  ±  1 1.4%  for  the  root  vector  difference.  ITiere  is  no  clear 
trend  with  object  type  or  orientation.  The  sources  of  the  estimation  bias  will  be  discussed 
in  the  next  Chapter. 
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Figure  4. 1 1 

Estimates  of  M\  component  of  the  dipole  moment  as  a  function  of 
Spheroid  depth  for  spiieroid  F  at  orientation  0  0.  Dots  are  estimates 
for  experimental  magnetic  field  maps.  Stars  are  for  theoretically  derived 
maps  (discussed  in  next  Chapter).  The  dashed  line  represents  the  dipole 
moment  component  calculated  for  the  object  type  and  orientation,  using 
the  model  of  Section  2.4.  Uncertainties  are  shown  where  available  and 
are  in  some  cases  smaller  than  the  data  points. 
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Figure  4.12 

Estimates  of  Mz  component  of  the  dipole  moment  as  a  function  of 
spheroid  depth  for  object  F  at  orientation  0  0.  Dots  ttre  estimates  for 
experimental  magnetic  field  maps.  Stars  are  for  theoretically  derived 
maps  (discussed  in  next  Chapter),  'llic  dashed  line  represents  the  dipole 
moment  component  calculated  for  the  object  type  and  orientation,  using 
the  model  of  Section  2,4,  Uncertainties  arc  shown  where  available  and 
are  in  some  cases  smaller  than  the  data  points. 
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Figure  4.13 

Estimates  of  A/3  component  of  the  dipole  moment  as  a  function  of 
spheroid  depth  for  object  F  at  orientation  0  0.  Dots  are  estimates  for 
experimental  magnetic  field  maps.  Stars  are  for  theoretically  derived 
maps  (discussed  in  next  Chapter).  The  dashed  line  represents  the  dipole 
moment  component  calculated  for  the  object  type  and  orientation,  using 
the  model  of  Section  2.4.  Uncertainties  are  shown  where  available  and 
are  in  some  cases  smaller  than  the  data  points. 
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Figure  4.14 

Estimates  of  Mi  component  of  the  dipole  moment  as  a  function  of 
spheroid  depth  for  object  F  at  orientation  90  175.  Dots  are  estimates 
for  experimental  magnetic  field  maps.  Stars  are  for  theoretically  derived 
maps  (discussed  in  next  Chapter).  The  dashed  line  represents  the  dipole 
moment  component  calculated  for  the  object  type  and  orientation,  using 
the  model  of  Section  2,4.  Uncertainties  arc  shown  where  available  and 
ai’e  in  some  cases  smaller  than  the  data  points. 
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Figure  4. 1 5 

Estimates  of  Mi  component  of  the  dipole  moment  as  a  function  of 
spheroid  depth  for  object  F  at  orientation  90  175.  Dots  are  estimates 
for  experimental  magnetic  field  maps.  Stars  are  for  theoretically  derived 
maps  (discussed  in  next  Chapter).  The  dashed  line  represents  the  dipole 
moment  component  calculated  for  the  object  type  and  orientation,  using 
the  model  of  Section  2.4.  Uncertainties  are  shown  where  available  and 
are  in  some  cases  smaller  than  the  data  points. 
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Figure  4.16 

Estimates  of  Afj  component  of  the  dipoie  moment  as  a  function  of 
spheroid  depth  for  object  F  at  orientation  90  175.  Dots  are  estimates 
for  experimental  magisetic  field  maps.  Stars  are  for  thccrctically  derived 
maps  (discussed  in  next  Chapter).  The  dashed  line  represents  the  dipole 
moment  component  calculated  for  the  object  type  and  orientation,  using 
the  model  of  Section  2.4,  Uncertainties  are  shown  where  available  and 
are  in  some  cases  smaller  than  the  data  points. 
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Figure  4.17 

Estimates  of  Mi  component  of  the  dipole  moment  as  a  function  of 
spheroid  depth  for  object  M  at  orientation  90  265.  Dots  are  estimates 
for  experimental  magnetic  field  maps.  Stars  are  for  theoretically  derived 
maps  (discussed  in  next  Chapter).  Tlie  dashed  line  represents  the  dipole 
moment  component  calculated  for  the  object  type  and  orientation,  using 
the  mode!  of  vSection  2.4.  Uncertainties  are  shown  where  available  and 
are  in  some  cases  smaller  than  the  data  points. 
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Figure  4.18 

Estimates  of  Mi  component  of  the  dipole  moment  as  a  ^'uiiction  of 
spheroid  depth  for  object  M  at  orientation  90  265.  Dots  are  estimates 
for  experimental  ma^  letic  field  maps.  Siisrs  are  for  theoretically  derived 
maps  (discussed  in  next  Chapter).  Tlic  d;'.:;h{  d  line  represents  the  dipole 
moment  component  calculated  for  the  object  type  and  orientation,  using 
the  model  of  Section  2.4.  Uncertainties  are  shown  where  available  and 
are  in  some  cases  smaller  than  the  data  points. 
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Figure  4.19 

Estimates  of  A/3  component  of  the  dipole  moment  as  a  function  of 
spheroid  depth  for  object  M  at  oricntatioii  90  265,  Dots  are  cstiinatcs 
for  eripeii.-nental  magnetic  field  maps.  Stars  are  for  theoretically  derived 
map.*;  (discussed  in  next  Chapter).  The  dashed  line  represents  the  dipole 
momc:n  ccunponent  calculated  for  the  object  type  and  orientation,  using 
the  model  of  Section  2.4.  Uncertainties  are  shown  where  available  and 
are  in  some  cases  smaller  than  the  data  points. 
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Figure  4.20 

Difference  between  Mi  component  estimate  and  theoretical  Mi,  as  a 
percentage  of  the  magnitude  of  theoretical  M  for  two  different  objects 
and  a  number  of  different  orientations.  Points  arc  labelled  with  a  letter 
for  object  type  and  two  numbers.  The  first  gives  the  polar  angle  and  the 
second  gives  the  azimuthal  angle  in  degrees.  Points  with  error  bars  arc 
unweighted  averages  of  a  number  of  estimates  at  different  depths.  Error 
bars  are  the  standard  deviations. 
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Figure  4.21 

Difference  between  Mi  component  estimate  and  theoretical  Mj,  as  a 
percentage  of  the  magnitude  of  theoretical  M  for  two  different  objects 
and  a  number  of  different  orientations.  Points  arc  labelled  with  a  letter 
for  object  type  and  two  numbers.  The  first  gives  the  polar  angle  and  the 
second  gives  the  azimuthal  angle  in  degrees.  Points  with  error  bars  are 
unweighted  averages  of  a  number  of  estimates  at  different  depths.  Error 
bars  arc  the  standard  deviations. 
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Figure  4.22 

Difference  between  component  estimate  and  theoretical  Mj,  as  a 
percentage  of  the  magnitude  of  theoretical  M  for  two  different  objects 
and  a  ntimber  of  different  orientations.  Points  are  labelled  with  a  letter 
for  object  type  and  two  number.  Tlie  first  gives  the  polar  angle  and  the 
second  gives  the  azimuthal  ang; '  in  degrees.  Points  with  error  bars  aie 
unweighted  averages  of  a  number  i)f  estimates  at  dilTerent  depths.  Error 
bars  ate  the  standard  deviations. 
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Figure  4.23 

Magnitude  of  the  vector  difference  between  estimated  M  and  theoretical 
M,  as  a  percentage  of  the  magnitude  of  theoretical  M  for  two  different 
objects  and  a  number  of  different  orientations.  Points  are  labelled  with  a 
letter  for  object  type  and  two  numbers.  The  first  gives  the  polar  angle  and 
the  second  gives  the  azimuthal  angle  in  degrees.  Points  with  error  bars 
are  unweighted  averages  of  a  number  of  estimates  at  different  depths. 
Error  bars  arc  the  standard  deviations. 
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5.  Performance  of  Location  Estimation 


5.1  Sources  of  Error  in  Location  Estimation 

In  the  previous  Chapter,  we  saw  that  the  location  algorithm  produced  estimates  of  both 
location  and  moment  components  that  deviated  from  the  measured  or  theoretical  values 
by  amounts  that  were  on  average  significantly  greater  than  the  experimental  precision 
(Figs.  4,20  to  4.23).  In  this  Chapter  we  will  analyse  the  sources  of  these  deviations  in  an 
attempt  to  identify  ways  to  minimize  or  eliminate  them. 

There  are  a  number  of  possible  sources  of  error  in  estimating  the  dipole  location  and 
moment  of  a  compact  object  when  using  the  algorithm  of  Section  2.2  with  a  total  field 
magnetometer.  These  may  be  grouped  in  the  following  manner: 

1.  Experimental  Errors  and  Uncertainties 

(a)  Sensor  Positional  Uncertainty  -  Deviations  of  the  sensor  from  its  assumed 
position  will  introduce  errors,  'fhe  dominant  error  term  is  due  to  deviations 
of  the  sensor  from  a  known  plane.  These  have  been  quantified  in  [10].  To 
summarize  briefly,  the  accuracy  of  dipole  parameter  estimates  was  found  to  be 
independent  of  tne  uncertainty  in  sensor  position  parallel  to  the  plane,  if  the 
latter  was  <5%  of  the  depth.  The  percentage  errors  in  both  position  and  dipole 
moment  component  estimates  were  roughly  equal  to  the  percentage  uncertainty 
(relative  to  the  depth)  in  sensor  position  orthogonal  to  the  plane,  if  the  latter  was 
<  5%  of  the  depth. 

(b)  Finite  Volume  of  Magnetic  Sensor  Head  -  ITie  algoiithm  assumes  that  a  magnetic 
measurement  is  taken  at  an  infinitesimal  point  in  space,  but  the  magnetic  sensor 
head  is  a  finite  cylinder  with  a  length  cT  ~  17  cm  and  a  diameter  of  ~  6  cm. 
The  measurement  point  has  been  assumed  to  be  the  geometric  center  of  the 
active  voIum<-.  If  the  field  gradient  is  small,  the  measured  field  will  be  an 
accurate  representation  of  the  true  field  at  that  point.  If  the  gradient  is  large, 
the  effective  measurement  point  will  be  shifted  from  the  geometric  center  by 
an  amount  which  is  a  function  of  the  field  gradient.  This  would  produce  larger 
errors  for  objects  and  orientations  whose  magnetic  field  maps  have  larger  field 
gradients.  The  effect  of  finite  sensor  head  volume  is  to  introduce  a  position 
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dependent  positional  uncertainty  in  all  three  directions.  Properly  quantifying 
the  uncertainty  would  require  extensive  modelling  of  the  operation  of  the  sensor 
head  in  relation  to  the  gradients  of  the  magnetic  field.s  of  typical  spheroids. 

(c)  Magnetic  Field  Measurement  Uncertainty  -  According  to  [10],  assumed  zero 
mean  additive  Gaussian  eirors,  the  accuracy  of  dipole  parameter  estimates  is  in¬ 
dependent  of  SNR  of  the  magnetic  field,  if  the  SNR  of  the  field  maximum  in  the 
measurement  plane  is  greater  than  about  25.  The  smallest  maximum  encoun¬ 
tered  in  this  study  was  ~  100  nl’.  The  quantization  error  for  the  magnetometer 
system  is  1  nT,  so  that  the  estimates  should  be  independent  of  SNR.  In  addition, 
the  spheroid  field  in  the  plane  can  have  a  large  gradient,  which  combined  with 
positional  uncertainty,  will  produce  errors  that  arc  neither  Gaussian  nor  zero 
mean.  Although  gradients  of  ~  2  nT/min  for  a  map  with  a  1000  nT  maximum 
field  and  ~  0.3  nT/mm  for  a  map  with  a  150  nT  maximum  field  were  typically 
encountered,  these  occurred  far  from  the  extrema.  Gradients  near  the  extrema 
arc,  as  one  would  expect,  small.  (This  is,  in  fact,  a  strength  of  the  extremum 
location  estimator).  Thus,  errors  due  to  magnetic  field  measurement  uncertainty 
arc  expected  to  be  negligible. 


2.  Location  Algorithm  Approximations 


(a)  Locating  the  Maximum  and  Minimum  of  the  Magnetic  Field  in  the  Measurement 
Plane  -  A  method  that  was  successful  with  theoretical  purely  dipolar  magnetic 
data  was  to  search  a  coarse  grid  (~  25  cm  spacing)  to  quickly  approximate  the 
maximum  and  minimum  positions  and  then  to  search  a  fine  grid  search  (~  2.5 
cm  spacing)  to  accurately  locate  the  extrema  [10].  For  this  study,  practical 
experimental  limitations  made  it  necessary  to  use  only  a  single  search  over  a 
fixed  grid  which  was  fine  in  the  A'l  dii'ection  (~  2  -  3  cm  spacing)  but  coarse  in 
the  Xj  direction  (~  5  -  10  cm  spacing),  'fhis  may  cause  an  error  which  is  depth 
dependent,  since  the  spatial  extent  of  the  field  is  a  function  of  depth  whereas 
the  grid  spacing  is  fixed.  In  addition,  the  position  of  the  extrema  are  located  by 
assuming  that  the  centroid  calculated  in  the  vicinity  of  the  extremum  coincides 
with  the  extremum  position.  For  skewed  peaks,  this  will  produce  an  error. 


(b)  Use  of  Total  Magnetic  Field  Measurements  -  ITe  algorithm  is  strictly  speaking 
correct  only  for  vertical  component  magnetic  field  measurements.  Total  mag¬ 
netic  field  values  approximate,  but  are  not  identical  to,  vertical  component  values 
in  geographical  regions  with  high  magnetic  inclinations,  such  as  the  location  of 
our  laboratory.  This  will  cause  a  bias  in  the  estimation  of  dipole  parameters,  but 
will  not  cause  an  increase  in  the  standard  deviation  of  the  estimates. 


(c)  errors  related  to  computer  implementation.  -  These  include  approximations  to 
inverse  functions  unique  to  the  algorithm  and  are  negligible  compared  to  other 
sources  of  error. 


3.  Deviation  of  Field  of  a  Compact  Object  From  Pure  Induced  Dipole 


UNCLASSIFIED 


DRES-SR-582 


UNCLASSIFIED 


61 


(a)  Higher  Order  Moments  -  Ferrous  objects,  when  placed  in  a  magnetic  field, 
such  as  that  of  the  earth,  acquire  induced  magnetization  whose  magnitude  and 
direction  is  a  function  of  the  primary  field  strength,  ferromagnetic  susceptibility 
of  the  body,  its  shape  and  its  orientation  with  respect  to  the  primary  field. 
This  produces  a  dipolar  field  that  can  be  attributed  to  a  dipole  moment  source 
located  at  the  object  center.  However,  for  all  objects  except  a  sphere,  there 
are  also  field  components  due  to  higher  order  moments.  ITiese  field  terms  are 
generally  smaller  than  and  decrease  with  distance  from  the  object  faster  than 
the  dipole  term.  At  close  distance,  the  fields  due  to  higher  order  moments  will 
add  a  distortion  term  to  the  dipole  field  and  lead  to  errors  in  dipole  parameter 
estimation.  For  the  algorithm  used  in  this  report,  the  effect  of  higher  multipole 
fields  is  to  shift  the  extrema  positions  and  alter  their  magnitude.  For  spheroids 
similar  in  shape  and  size  to  artillery  shells  at  close  range  in  geometries  similar 
to  those  of  these  experiments  (depth  4  times  the  maximum  dimension  and  18 
times  the  minimum  dimension),  the  maximum  octupole  field  varies  from  3  to 
8%  of  the  maximum  dipole  field  [8]. 

(b)  Remnant  Magnetiz.ation  -  Ferrous  objects  may  also  have  remnant  magnetization 
which  is  a  function  of  the  metallurgical  properties  of  the  object  as  well  as  its 
thermal,  mechanical  and  geomagnetic  history.  Remnant  magnetization  is  very 
difficult  to  model  since  its  strength  and  direction  are  often  unknown  for  an 
individual  object  and  may  vary  from  object  to  object  of  the  same  shape  and 
size.  Remnant  magnetization  will  in  general  not  be  purely  dipolar.  The  dipole 
component  of  the  remnant  magnetization  will  add  vectorally  to  the  induced 
moment,  which  may  make  it  difficult  to  use  the  estimated  total  dipole  moment  to 
classify  the  object.  Higher  order  vnultipoles  will  cause  a  distortion  to  the  dipole 
field,  just  as  those  from  the  induced  magnetization  do,  which  will  increase 
the  error  in  estimating  the  source  location  and  dipole  moment.  Fortunately, 
past  evidence  has  suggested  that  remnant  magnetization  can  be  neglected  for 
some  applications,  including  ordnance  detection  |7J.  Even  when  present  in 
artillery  shel's,  some  evidence  indicates  that  remnant  magnetization  is  small 
and  is  similar  for  individual  shells  of  the  same  type  [8J. 

Of  the  non-negligible  errors  mentioned  above,  positional  uncertainties  will  random  in 
nature  and  will  contribute  to  the  precision  of  the  estimates  but  not  the  accuracy  (bias).  The 
maximum  uncertainties  in  position  were  roughly  1.4  cm  for  Xi  (determined  mainly  by 
the  finite  measureinenl  time  of  the  magnetometer),  1  cm  in  A'z  (dictated  by  swerving  of 
the  magnetometer  cart)  and  0.2  to  0.5  cm  in  .V3  (caused  by  fluctuations  in  flatness  of  the 
measurement  table).  Since  the  uncertainty  in  both  Xi  and  X2  was  <  5%  of  the  depth  for 
all  depths,  [10]  implies  that  these  uncertainties  should  have  a  negligible  contribution  to  the 
overall  precision.  To  compare  the  measured  precision  with  the  precision  attributable  to  the 
positional  uncertainty,  we  can  average  the  component  preci.sion  of  the  location  and  moment 
components  of  Tables  VII  and  VIIl  to  produce  Table  X.  It  is  seen  that  the  precision  of 
estimated  location  components  is  in  reasonable  agreement  with  [10],  that  is,  the  average 
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precision  in  the  components  is  roughly  equal  to  the  percentage  uncertainty  in  X3.  The 
precision  of  the  estimated  moment  components  is  slighly  higher  than  the  uncertainty  in 
X'i  but  is  still  in  reasonable  agreement  with  [10].  In  any  case,  it  is  clear  that  a  substantial 
fraction,  if  not  the  major  portion,  of  the  uncertainty  (precision)  in  dipole  location  and 
moment  vector  estimates  is  due  to  positional  uncertainty. 

ITie  bias  in  estinnation  in  the  dipole  location  and  moment  vector  estimates  must  then  be 
due  to  the  remaining  five  non-negligible  sources  of  error.  We  shall  next  study  the  error  in 
estimation  using  magnetic  field  maps  generated  by  computer  using  a  mathematical  model 
to  resolve  some  of  the  remaining  sources. 


5.2  Location  Estimates  for  Theoretical  Total  P'ield  Data 

If  the  location  algorithm  is  applied  to  magnetic  field  data  generated  by  the  mathematical 
model  of  Section  2.4,  errors  due  to  positional  and  magnetic  measurement  uncertainties 
are  eliminated.  Also,  the  true  magnetic  dipole  moment  is  available  to  compare  with  the 
estimated  dipole  moment,  s' nee  there  is  no  remnant  magnetization.  It  will  then  be  possible 
to  investigate  bias  in  estimation  due  to  finite  volume  of  the  sensor  head,  finding  the  field 
maximum  and  minimum  in  the  plane,  higher  order  moments  and  the  use  of  total  field 
measurements. 

Examples  of  typical  magnetic  field  maps  obtained  using  the  model  are  shown  in  Fig¬ 
ures  5.1  to  5.4.  These  should  be  compared  with  Figures  4.1  to  4.4.  The  shapes  and  peak 
field  values  of  the  experimental  and  theoretical  maps  are  very  similar  as  one  would  ex¬ 
pect  given  the  similarity  between  theoretical  and  experimental  dipole  moments  (Figs.  4.20 
to  4.22). 

The  location  estimation  algorithm  was  applied  to  the  same  object  types  and  orientations 
at  depths  similar  to  those  that  were  used  in  the  exfjerimental  measurements.  The  results 
are  shown  in  Tables  XI  and  XII.  These  should  be  compaied  with  Tables  V  and  VI.  Also, 
in  Figs.  4.11  to  4.19,  estimated  dipole  moment  components  for  magnetic  total  field  maps 
generated  by  computer  using  the  mathematical  model  are  shown  in  comparison  to  those 
estimated  from  experimental  field  maps. 

The  rms  error,  averaged  (unweighted)  over  all  cases  in  Tables  V,  XI,  VI  and  XII,  and  its 
standard  deviation  are  used  as  a  rough  measure  of  estimation  pciformance  for  purposes  of 
comparison.  From  the  tables,  the  average  enor  in  estimation  of  the  location  vector  is  0.04 
m  ±  0.01  m  for  computer  generated  data  versus  0.05  m  ±  0.02  m  for  experimental  data. 
Tlie  average  error  in  estimation  of  the  dipole  moment  vector  is  22  %  ±  6  %  (expressed  as  a 
percentage  of  the  moment  magnitude)  for  computer  generated  data  versus  27  %  ±  10  %  for 
experimental  data.  Estimation  errors  for  both  experimental  and  computer  generated  data 
are  significantly  larger  than  the  average  precision,  which  from  Tables  VII  and  VTII  is  0.003 
m  ±  0.002  m  for  the  location  vector  and  4.1  ±  3.0  %  for  the  moment  vector.  Tlie  error  in 
estimation  of  the  location  vector  is  largely  independent  of  whether  the  data  is  experimental 
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Figure  5.1 

"Mieorcticai  total  magnetostatic  held  map  versus  position  in  a  plane  for 
spheroid  F  at  0  -  0°.  Spheroid  parameters  arc  shown  at  right.  Ambient 
(earth’s)  field  magnitude  has  been  subtracted. 
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Figure  5.2 


Tncorerical  (otai  magnccosiaiic  neid  map  versus  pc  (ion  in  a  pianc  lor 
spheroid  F  at  0  =  30°  and  <i>  ~  175°.  Spheroid  parameters  arc  shown  at 
right.  Ambient  (earth’s)  held  magnitude  has  been  subtracted. 
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Figure  5.3 

Theoretical  total  magnetostatic  field  map  versus  position  in  a  plane  for 
spheroid  F  at  ^  -  60®  and  <!>  =  1 75®. Spheroid  parameters  arc  shown  at 
right.  Ambient  (earth’s)  field  magnitude  has  been  subtracted. 
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Figure  5.4 

Theoretical  total  magnetostatic  field  map  versus  position  in  a  plane  for 
spheroid  F  at  0  =  90°  and  ^  =  175°.  Spheroid  parameters  arc  shown  at 
right.  Ambient  (earth’s)  field  magnitude  has  been  subtracted. 
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or  theoretical.  Trends  in  overestimation  or  underestimation  of  the  depth  are  similar  to  those 
discussed  in  Chapter  4  for  the  experimental  data.  Dipole  moment  estimation  errors  for 
computer  generated  data  are  slightly  smaller  than  those  for  experimental  data.  This  implies 
that  the  positional  error  in  these  experiments  has  little  influence  on  location  and  moment 
vector  estimation.  This  alsi  suggests  that  remnant  magnetization  in  the  spheroids  must 
be  small  compared  to  the  induced  magnetization,  since  the  theoretical  moments  have  been 
used  for  comparison  to  the  estimates  for  both  the  theoretical  and  experimental  magnetic 
maps. 

To  further  investigate  trends  as  a  function  of  object  shape,  orientation  and  depth,  the 
rms  estimation  error  for  vertical  and  horizontal  type  M  and  F  spheroids  were  plotted  for 
a  number  of  depths.  In  addition,  the  rms  estimation  error  for  a  spheroid  that  was  neither 
vertical  nor  horizontal  and  a  sphere  were  also  plotted.  These  are  shown  in  Figs.  5.5  and 
Figs.  5.6. 

The  figures  reveal  a  number  of  {Joints.  The  error  in  estin.,.tion  of  the  location  vectors 
is  generally  bigger  than  the  precision  in  estimation  (Fig.  4.5)  for  both  computer  generated 
and  experimental  data.  Since  the  bias  for  both  data  sets  is  similar,  this  suggests  that  the 
finite  volume  of  the  magnetometer  sensor  head  is  not  a  significant  source  of  error.  (Because 
the  field  gradient  decreases  with  distance  from  the  source,  the  tight  geometry  of  these 
experiments  is  probably  a  worst  case  for  enors  due  to  finite  sensor  volume.)  Generally 
speaking,  the  eaor  shows  some  variation  with  depth  but  it  appears  almost  random.  'Ihe 
error  also  appears  to  be  independent  of  object  shap)e.  The  error  appears  to  depend  slightly 
on  the  object  orientation  and  the  same  orientation  yields  similar  errors  for  different  object 
shapes.  Specifically,  the  90  90  orientation  yields  the  smallest  erroi  at  low  depths  (<  0.7m) 
and  the  10  0  orientation  is  smallest  at  greater  depths.  Tlie  90  175  orientation  yields  the 
largest  error  for  all  depths.  Since  the  error  for  the  sphere  is  in  the  middle  of  the  range 
of  errors,  higher  order  moments  do  not  seem  to  have  a  large  effect  on  the  estimation  of 
location.  The  variation  of  ei  ror  with  depth  suggests  that  the  main  conuibutor  to  the  error  is 
estimation  of  the  position  of  the  field  extrema  and  the  use  of  total  field  in  place  of  vertical 
component  magnetic  data,  since  both  are  de{Jendent  on  depth  for  a  fixed  sampling  grid. 
These  two  sources  of  error  are  further  indicated  because,  as  has  been  seen,  both  errors 
should  show  a  de|)cndence  on  the  object  orientation. 

The  spread  of  the  moiiieni  eslimalion  error  curves  is  comparable  to  precision  of  moment 
estimates  (Fig  4.6),  but  there  is  a  clear  bias  in  estimation  for  all  curves.  The  error  shows 
a  small,  seemingly  random  dependence  on  depth  except  for  the  F  90  175  error  at  depths 
<  0.5  m  which  increases  with  decreasing  depth.  The  error  for  the  sphere  is  smaller  than 
any  other  object  or  orientation  except  for  F  90  175  at  depths  >  0.6  m.  This  suggests  that 
the  presence  of  higher  order  multipole  moments  is  significant  in  di{X)le  moment  estimation, 
but  the  effect  is  small  compared  to  the  overall  bias  for  the  error  curves.  Since  the  shapes 
of  magnetic  field  maps  are  a  function  of  object  shape  and  orientation,  it  is  not  possible  to 
separate  the  effects  of  higher  order  mullipoles  from  the  error  in  extrema  {josilion  estimation. 
The  imjxjrtant  point  for  object  classification  is  that  the  momeni:  estimates  for  a  given  object 
and  orientation  vary  by  between  2  and  4%  of  the  moment  magnitude  as  the  depth  varies. 
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Fi  cure  5-5 

rms  location  vector  estimation  error  for  computer  generated  magnetic 
total  field  maps  for  two  spheroids  at  various  oiicntations  and  a  sphere. 
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Figure  5,6 

rms  moment  vector  estimation  error  for  computer  generated  magnetic 
total  ,.eld  maps  for  two  spheroids  at  various  orientations  and  a  sphere. 
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The  bias  is  unimportant,  since  the  estimated  moments  can  be  used  in  the  design  set.  This 
suggests  that  the  limiting  misclassification  error  will  be  between  roughly  1 .4  and  2.4%  [11] 
if  the  experimental  error  is  substantially  less  than  the  algorithmic  and  approximation  en'ors. 
For  the  present  experimental  precision  (4.1  ±  3.0%  ),  the  misclassification  error  should  be 
roughly  between  2  and  5%  (assuming  quadrature  error  summation)  if  the  design  set  were 
experimentally  obtained. 

For  both  the  location  and  dipole  moment  estimation  the  F  90  175  curves  do  not  follow 
the  same  mends  as  the  others.  T!ie  location  enor  is  higher  tb^.n  the  other  errors  while  the 
moment  curve  has  a  sharp  depth  dependency  not  seen  in  the  other  moment  errors.  The 
reason  for  this  is  not  known. 

It  is  initially  tempting  to  think  that  one  can  estimate  the  remnant  dipole  moment  by 
subtracting  the  dipole  moments  estimated  for  experimental  maps  from  those  estimated 
for  theoretical  maps  with  object  type,  orientation  and  depth  fixed  (Figs.  4.11  to  4.19). 
Averaging  over  depth  would  then  give  the  remnant  magnetization  for  a  given  object  and 
orientation  (which  should  be  related  to  other  orientations  by  a  simple  Euler  rotation  matrix). 
Unfortunately,  as  we  have  just  seen,  the  bias  in  estimation  varies  with  object  orientation. 
This  is  presumably  due  to  the  differing  shapes  of  the  magnetic  maps  which  have  different 
errors  when  estimating  the  extrema  positions.  This  variation  is  small  but  is  clearly  of  the 
same  order  of  magnitude  as  the  deviation  between  experimental  and  theoretical  estimates. 
A  simple  subtraction  of  experimental  and  theoretical  moment  estimates  would  not  take 
into  account  the  difference  in  bias  due  to  the  differing  moments  of  the  experimental  and 
theoretical  data.  There  may  be  ways  to  more  precisely  estimate  the  contribution  of  remnant 
magnetization  to  the  total  dipole  moment.  For  example,  one  might  use  theoretical  moments 
to  derive  the  error  bias  as  a  function  of  moment  component  estimate.s.  Tl.is  could,  in  turn, 
be  applied  as  a  correction  to  the  experimental  estimates.  Then,  one  could  subtract  the 
theoretical  dipole  moments  from  the  corrected  experimental  estimates.  It  is  hoped  that  this 
can  be  investigated  in  the  future. 


5.3  Location  Estimates  for  Theoretical  Vertical  Component  Field  Data 

From  the  previous  section,  it  was  seen  that,  if  we  exclude  experimental  errors  (positional 
uncertainty  and  remnant  magnetization),  there  are  three  other  possible  sources  of  eiTor  that 
can  contribute  significantly  to  the  estimation  error,  lli  sc  are  the  estimation  of  extrema 
positions,  the  use  of  toial  field  magnetic  data  and  the  presence  of  higher  order  moments. 
To  eliminate  the  effect  of  using  the  total  magnetic  field  of  the  spheroid  as  ati  approximation 
to  the  vertical  component  of  the  magnetic  field,  the  location  algorithm  'vas  applied  to 
magnetic  field  maps  of  the  vertical  magnetic  component  Uh)-  The  rms  estimation  errt 
was  calculated  for  the  same  objects,  orientations  and  depths  as  in  Pigs.  5.5  and  5.6.  'iliese 
are  shown  in  Figs.  5.7  and  Figs.  5.8. 

ITie  error  in  estimation  of  the  location  vector  is  generally  slightly  smaller  for  a  given 
case  than  when  total  field  data  are  used.  Other  trends  are  very  .similar  to  the  total  field  data. 
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Figure  5.7 


rms  location  vector  estimation  error  for  computer  generated  magnetic 
vertical  component  field  maps  for  two  spheroids  at  various  orientations 
and  a  sphere. 
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Figure  5.8 

rms  moment  vectc  estimation  error  for  computer  generated  magnetic 
vertical  component  field  maps  for  two  spheroids  at  various  orientations 
and  a  sphere. 
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Again  the  depth  dependency  is  small  and  almost  random  in  nature.  The  error  for  the  sphere 
is,  as  before,  not  the  lowest  for  all  cases.  The  90  175  case  again  yields  the  highest  error  at 
all  depths,  except  for  M  90  90  at  0.5  m.  It  is  not  clear  why  the  error  for  this  latter  case  is 
anomalously  high. 

The  spread  in  moment  cst  mation  error  for  vertical  component  data  is  similar  to  that  of 
the  total  field  data,  but  the  bias  for  the  former  is  significantly  less.  This  suggests  that  the 
use  of  total  field  data  introduces  a  bias  in  estimation  for  e  dipole  moment.  As  for  the  total 
field  data,  the  error  for  vertical  component  data  is  roughly  constant  with  depth  for  a  given 
object  and  orientation.  The  error  in  estimation  of  the  sphere  is  significantly  less  than  all 
other  cases  except  M  90  90  and  M  0  0  at  0.6  m  depth,  which  indictates  that  higher  order 
moments  still  contribute  somewhat  to  the  error.  To  crudely  estimate  the  contribution  of 
the  various  sources  of  error,  we  note  that  the  bias  in  estimation  of  the  sphere  for  vertical 
component  data  is  roughly  2  to  6%  of  the  moment  magnitude,  whereas  for  total  field  data 
it  is  20%.  This  suggests  that  the  error  due  to  the  total  field  approximation  varies  from 
14  to  18%  of  the  moment  magnitude.  Error  due  to  estimation  of  extrema  positions  and 
higher  order  moments  then  varies  from  0  to  8%  of  the  moment  magnitude  which  is  in  rough 
agreement  with  that  for  the  total  field  data.  This  seems  reasonable,  since  in  Section  5.1  it 
was  noted  the  maximum  octupole  field  varies  from  3  to  8%  of  the  maximum  dipole  field 
for  spheroids  in  geometries  similar  to  ours. 
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Depth 

Range  of  Uncertainty 

Average  Location 

Average  Moment 

in  X3 

Component  Precision 

Component  Precision 

(m) 

(%  of  Depth) 

(%  of  Depth) 

(%  of  Moment  Magnitude) 

0.390 

0.5- 1.3 

2.0 

4.7 

0.438 

0.4- 1.1 

1.0 

2.5 

0.500 

0.4- 1.0 

1.9 

- 

0.525 

0.4- 1.0 

3.9 

- 

0.565 

0.4  -  0.9 

0.2 

0.3 

0.598 

0.3  -  0.8 

0.8 

2.4 

0.698 

0.3  -  0.7 

0.2 

0.7 

0.802 

0.2  -  0.6 

0.5 

1.8 

Table  X 

Comparison  of  uncertainty  in  Xj  measurement  with  precision  of  esti¬ 
mated  location  and  dipole  moment  vector  components.  Precision  of 
location  (moment)  components  is  the  average  of  the  precision  of  the  3 
location  (moment)  components  for  a  given  depth,  taken  from  Tables  VII 
and  VIII. 
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Spheroid 

IVpe, 

Orientation 

h 

Xoi 

(m) 

/leasuret 

A'o2 

(m) 

i 

Xo3 

(m) 

E 

Ao, 

(m) 

istimated 

Ao2 

(m) 

Ao3 

(m) 

FOO 

0.000 

0.000 

0440 

-0.004 

0.000 

0.399 

0.000 

0.000 

0.500 

-0.015 

0.000 

0.477 

0.000 

0.000 

0.600 

-0.007 

0.000 

0.568 

0.000 

^000 

0.802 

-0.014 

0.000 

0.772 

0.000 

0.000 

1.000 

-0.022 

0.000 

^0.978 

F90  175 

0.000 

0.000 

0.440 

-0.028 

0.004 

0.500 

0.000 

0.000 

0.500 

-0.031 

0.010 

0.546 

0.000 

0.000 

0.525 

-0.037 

0.010 

0.546 

0.000 

0.000 

0.600 

-0.037 

-0.003 

0.579 

0.000 

0.000 

0.655 

-0.031 

0.004 

0.673 

0.000 

0.000 

0.802 

-0.045 

0.003 

0.818 

0.000 

0.000 

1.000 

-0.051 

0.003 

1.018 

M  90  265 

0.000 

0.000 

0.440 

-0.015 

0.000 

0.450 

0.000 

0.000 

0.500 

-0.009 

0.000 

0.500 

0.000 

0.000 

0.600 

-0.011 

0.000 

0.583 

0.000 

0.000 

0.802 

-0.022 

0.000 

0.789 

0.000 

0.000 

1.000 

-0.019 

0.000 

0.987 

M  10  175 

0.000 

0.000 

0.500 

0.010 

0.000 

0.472 

F60  175 

0.000 

0.000 

0.525 

-0.041 

0.004 

0.546 

F30  175 

0.000 

0.000 

0.525 

-0.041 

-0.003 

0.521 

F  60  265 

0.000 

0.000 

0.525 

-0.040 

-0.021 

0.553 

F  30  265 

0.000 

0.000 

0.525 

-0.018 

-0.031 

0.504 

Table  XI 

Location  vector  estimation  results  using  magnetic  field  data  generated 
from  the  computer  model  of  Section  2.4. 
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ITieoretical 

Estimated 

spheroid 

depth 

Mx 

A/2 

A/3 

Mx 

Mi 

A/3 

type, 

orientation 

(m) 

(Am^) 

(A-m^) 

(Am^) 

(A-m^) 

(A-m^) 

(A.m^) 

FOO 

0.440 

0.0649 

0.0000 

1.0685 

0.2685 

0,0000 

0.9073 

0.500 

0.0649 

0.0000 

1.0685 

0.3072 

0.0000 

1.0144 

0.600 

0.0649 

0.0000 

1.0685 

0.2965 

0.0000 

0.9535 

0.802 

0.0649 

0.0000 

1.0685 

0.3081 

0.0000 

0.9593 

1.000 

0.0649 

0.0000 

^1.0685 

0.3219 

0.0000 

0.9881 

F  90  175 

0.440 

0.3247 

-0.0227 

0.2123 

0.4640 

-0.0361 

0.2326 

0.500 

0.3247 

-0.0227 

0.2123 

0.4214 

-0.0360 

0.2137 

0.525 

0.3247 

-0.0227 

0.2123 

0.4362 

-0.0470 

0.2208 

0.600 

0.3247 

-0.0227 

0.2123 

0.3707 

-0.0306 

0.1864 

0.655 

0.3247 

-0.0227 

0.2123 

0.3614 

-0.0280 

0.1828 

0.802 

0.3247 

-0.0227 

0.2123 

0.3600 

-0.0162 

0.1812 

1.000 

0.3247 

-0.0227 

0.2123 

0.3615 

-0.0231 

0.1820 

M  90  265 

0.440 

0.0579 

0.0000 

0.1893 

0.1033 

0.0000 

0.1802 

0.500 

0.0579 

0.0000 

0.1893 

0.0963 

0.0000 

0.1685 

0.600 

0.0579 

0.0000 

0.1893 

0.0889 

0.0000 

0.1575 

0.802 

0.0579 

0.0000 

0.1893 

0.0923 

O.OCOO 

0.1655 

1.000 

0.0579 

0.0000 

0.1893 

0.0949 

0.0000 

0.1674 

M  10  175 

0.500 

0.1045 

-0.0041 

0.4548 

0.1106 

0.0000 

0.3585 

F  60  175 

0.525 

0.6291 

-0.0494 

0.5393 

0.7586 

-0.0649 

0.5364 

F  30  175 

0.525 

0.4992 

-0.0380 

0.9674 

0.6788 

-0.0199 

0.9528 

F  60  265 

0,525 

0.0987 

0.3864 

0.4363 

0.1423 

0.4281 

0.4861 

F  30  265 

0.525 

0.0977 

0.3750 

0.8644 

0.2140 

0.2951 

0.8032 

Table  XII 

Dipole  moment  vector  estimation  results  using  magnetic  field  data  gen¬ 
erated  from  the  computer  model  of  Section  2.4. 
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6.  Performance  of  Pattern  Classification 


Identification  of  objects  from  the  estimated  dipole  moments  was  achieved  using  the 
continuous  parameter  pattern  (CP)  classifier  of  Section  2.3.  Following  [11],  the  design  set 
foreach  object  consisted  ofadipolemomentforeach  15®  increment  in  0  and  c6,0  <  0  <  90°, 
0<<j><  360°. 

T1.C  design  set  for  one  object  requires  225  dipole  moment  vectors.  At  8  -  15  minutes 
per  dipole  moment  measurement  for  the  present  instrument,  a  design  set  using  njornents 
estimated  from  experimental  data  would  require  approximately  30  to  55  measurement  hours 
per  object.  Since  about  8  to  10  objects  are  needed  to  reasonably  represent  the  range  of 
typical  unexploded  ordnance  that  is  available,  between  240  and  550  hours  of  measurement 
time  would  be  required  in  total.  On  the  other  hand,  a  dipole  moment  design  set  derived 
from  Equation  2.33  of  the  model  of  Section  2.4  could  be  obtained  in  a  few  minutes.  Since 
our  aim  was  to  test  the  proof  of  concept  of  real-time  location  and  identification  of  objects, 
in  the  interests  of  expediency  we  chose  to  use  the  theoretical  design  set.  In  doing  so,  it  was 
recognized  that  a  systematic  difference  between  estimated  and  theoretical  dipole  moments, 
as  we  have  seen  in  the  previous  Section,  would  increase  the  misclassification  error  rate  over 
what  would  be  achieved  with  an  experimentally  derived  design  set. 

The  design  set  used  in  the  classification  is  shown  in  Table  XlII. 

Results  of  the  classification  using  the  experimentally  estimated  moments  are  given  in 
Table  XIV  for  spheroid  F  and  in  Table  XV  for  spheroid  M.  For  each  case,  the  classification 
algorithm  provided  a  ranked  list  of  all  classes  and  the  closest  distance  from  each  class  to 
the  test  vector.  Here,  only  the  highest  ranked  (most  likely  class)  has  been  shov  n  for  each 

OOCIA 

SpheroidF  wasquite  successfully  classified  with  a  misclassification  rate  of  3/29  (10.3%). 
The  orientation  angles  of  the  spheroid  were  also  estimated  quite  well.  For  0,  all  estimates 
for  correctly  classified  cases  were  within  10°  of  the  true  angle  and  roughly  46%  were  within 
5°.  For  (^,  all  estimates  for  correctly  classified  cases  were  within  15°  of  the  true  angle, 
roughly  92%  were  within  10°  and  69%  were  within  5°. 

At  first  glance,  results  were  less  encouraging  for  Spheroid  M.  None  of  the  cases  were 
classified  correctly.  However,  the  spheroid  was  consistently  classified  as  class  number  3. 
All  estimates  of  0  for  the  90  265  orientation  were  within  21°  of  each  other  and  5/6  were 
within  5°  of  each  other.  4/6  estimates  of  <f>  for  the  90  265  orientation,  all  at  the  same  depth. 
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were  svithin  5°  of  each  other.  The  remaining  2/6  were  not  even  close  to  one  another  or  the 
others. 

Further  light  can  be  shed  on  the  classification  error  by  examining  the  individual  mis- 
classifications  in  Table  XVI.  Here  we  have  listed  classification  results  for  each  case  ranked 
in  descending  order  down  to  and  including  those  for  the  true  object.  The  distance  measure 
presented  in  column  7  is  (Equation  2.16). 

Of  the  misclassifications  of  object  6,  the  correct  object  is  ranked  second  4/7  times.  The 
estimates  of  angles  have  a  larger  spread  for  object  3  than  for  the  correct  classifications  of 
object  4.  'rhis  suggests  that,  although  object  3  is  chosen  as  the  most  likely  class  when 
experimental  data  for  object  6  is  presented,  the  computer  generated  dipole  moments  of 
object  3  are  not  as  good  an  approximation  to  the  experimental  moments  of  object  6  as 
computer  generated  moments  of  object  4  are  to  experimental  moments  of  object  4.  The 
consistency  with  which  object  6  is  classified  as  object  3  can  be  explained  by  noting  that 
class  3  and  class  6  objects  are  the  closest  in  shape  of  the  objects  in  the  design  set.  A  class  3 
spheroid  has  a  length  of  22.5  cut  and  a  width  of  9  cm,  while  a  class  6  spheroid  has  a  length 
of  23.6  cm  and  a  width  of  11.8  cm.  A  small  amount  of  remnant  magnetization  together 
with  the  previously  described  errors  in  estimating  the  dipole  moment,  can  easily  account 
for  the  misclassification. 

Object  4  is  misclassified  3  times  as  a  sphere.  In  addition,  object  6  is  ranked  ahead 
of  it  in  all  3  cases.  Since  object  4  is  classified  correctly  the  remainder  of  the  time,  these 
misclassifications  are  likely  a  reflection  of  the  error  in  estimating  the  dipole  moment  in 
the  depth  range  between  0.5  and  0.6  m.  Tliat  rant‘;e  exhibits  the  largest  deviation  between 
experimental  and  computer  generated  moments  for  Spheroid  F  at  0  0  (Figs.  4.1 1  to  4,13). 

The  consistency  or  repeatability  of  classification  is  really  more  important  than  the  correct 
choice  of  object,  since  if  object  “a”  is  always  classified  as  object  “b”,  one  can  merely  switch 
the  as.signments  upon  completing  the  classification.  (This  assumes,  of  course,  that  object 
“b”  is  always  classified  as  something  else!)  If  we  thus  consider  the  experimental  dipole 
moments  for  Spheroid  M  to  be  most  closely  represented  by  the  computer  generated  moments 
of  the  class  3  spheroid,  the  misclassification  rate  is  3/36  (8.3%). 

I  inally,  if  we  consider  only  the  two  objects  for  which  experimental  measurements  were 
made,  we  see  that  only  4/36  cases  (11.1%)  were  classified  incorrectly. 

Overall,  the  misclassification  rate  is  encouraging  l  ut  improvement  is  desirable.  Tlie 
misclassification  rate  of  the  pattern  classifier  is  integrally  linked  to  the  error  in  estimation  of 
the  dipole  moment  by  the  location  algorithm.  Tltere  are  a  number  of  things  that  can  be  done 
to  improve  the  misclassification  rate.  One  obvious  solution  is  to  use  a  design  set  based  on 
estimated  dipole  moments  from  experimental  data.  Although  this  would  be  fime  consuming 
using  the  instrument  of  this  investigation,  such  a  design  set  will  be  obtained.  Further,  a  new 
multisensor  fluxgate  gradiometer,  which  we  are  presently  developing,  will  allow  location 
and  dipole  moment  estimation  without  moving  the  sensor.  If  successful,  this  should  speed 
up  design  set  acquisition  by  at  least  a  factor  of  five.  Better  extrema  estimation  will  improve 
the  misclassification  rate.  This  might  be  achieved  by  locating  the  true  peaks  in  the  field  map 
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as  opposed  to  locating  the  cenuoid.  One  way  to  do  this  is  to  regularize  the  sampling  grid, 
so  that  two  dimensional  interpolation  routines  can  be  used  to  more  accurately  pinpoint  the 
extrema.  If  depth  dependence  in  the  error  bias  cannot  be  eliminated,  design  sets  could  be 
stored  for  various  depths  to  minimize  the  depth  dependent  estimation  bias.  Problems  with 
this  approach  include  the  need  for  substantial  amounts  of  microprocessor  memory  and  the 
time  to  collect  a  very  large  design  set.  A  model  for  bias  as  a  function  of  estimated  moment 
component  may  be  developed.  This  could  then  be  applied  as  a  correction  to  the  estimated 
moments  prior  to  classification.  These  items  will  be  the  subject  of  later  investigations. 
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unclassieied 


Class  Number 

Ubel 

rt(m) 

e 

T^pc  of 

Unexploded  Ordnance 
Resembled 

0 

0.01 

3.75 

20mm  shell 

1 

0.02 

3.4 

40mni  shell 

2 

0.03 

2.5 

60mm  shell 

3 

0.045 

2.5 

81mm  shell 

4 

Spheroid  F 

0.0525 

3.5 

105mm  shell 

5 

0.09 

3.5 

155mm  shell 

6 

Spheroid  M 

0.05896 

2.0 

7 

0.10 

1.0 

sphere 

Table  XIII 

Objects  used  in  design  set  for  classification  studies.  Spheroid  parameters 
■ire  defined  in  Section  2.4.  Labels  refer  to  the  objects  for  which  exper- 
mental  data  is  available.  Last  column  shows  the  unexploded  ordnance 
type  that  the  spheroid  is  intended  to  loosely  represent. 
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Class 

Number 

0 

n 

4> 

n 

depth 

(m) 

Estimated 

0 

O 

Estimated 

<t> 

n 

A 

0 

- 

4 

5.05 

193.08 

4 

0 

- 

4 

6.83 

192.27 

4 

0 

- 

4 

8.35 

193.35 

4 

0 

- 

4 

7.20 

189.43 

4 

0 

- 

4 

6.99 

207.65 

4 

0 

- 

4 

6.49 

192.72 

4 

0 

- 

4 

7.53 

190.53 

4 

0 

- 

4 

7.37 

193.12 

4 

0 

- 

0.598 

4 

6.87 

191.14 

4 

0 

- 

4 

8.90 

173.51 

4 

0 

- 

4 

10.31 

178.30 

4 

0 

- 

4 

5.91 

221.92 

4 

0 

- 

7 

0.00 

15.00 

4 

0 

- 

7 

0.00 

15.00 

4 

0 

- 

7 

15.00 

0.00 

4 

90 

175 

0.655 

4 

84.18 

0.00 

4 

90 

175 

4 

85.80 

174.17 

4 

90 

175 

0.44 

4 

89.29 

0.71 

4 

90 

175 

4 

85.25 

12.59 

4 

90 

175 

4 

87.82 

6.13 

4 

90 

175 

4 

86,52 

4.01 

4 

90 

175 

4 

85.79 

0.00 

4 

90 

175 

0.39 

4 

89.32 

169.79 

4 

90 

175 

iia 

4 

84.89 

0.00 

4 

90 

175 

119 

4 

85.98 

0.00 

4 

60 

175 

0.525 

4 

59.61 

176.44 

4 

175 

4 

39.37 

175.46 

4 

60 

265 

0.525 

4 

52.62 

259.39 

4 

30 

265 

4 

32.13 

256.86 

Table  XIV 

Results  of  classifying  dipole  moments  of  spheroid  F  estimated  from 
experimental  magnetic  field  maps  using  the  CP  classifer  and  the  design 
set  of  Table  XIIl. 
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Class 

Number 

0 

n 

<l> 

C) 

depth 

(m) 

Closest 

Class 

Number 

Estimated 

0 

n 

Estimated 

<!> 

n 

6 

90 

265 

0.58 

3 

70.36 

19.13 

6 

90 

265 

0.468 

3 

86.43 

254.58 

6 

90 

265 

0.438 

3 

68.56 

126.30 

6 

90 

265 

0.438 

3 

70.60 

124.40 

6 

90 

265 

0.438 

3 

71.40 

123.59 

6 

90 

265 

0.438 

3 

67.64 

127.36 

6 

10 

175 

0.50 

3 

3.09 

29.08 

Table  XV 

Results  of  classifying  dipole  moments  of  spheroid  M  estimated  from 
experimental  magnetic  field  maps  using  the  CP  classifer  and  the  design 
set  of  Table  XIII. 


UNCLASiilHED 


DRES-SR-582 


DRES-SR-582 


7.  Conclusions 


An  improved  “smart”  microprocessor  controlled  magnetometer  which  can  accurately 
locate  and  identify  contpact  ferrous  objects  in  real-time  has  been  described.  The  instrument, 
which  is  person-portable,  consists  of  a  cesium  vapour  magnetometer  mounted  on  a  cart  with 
a  wheel-mounted  optical  encoder,  a  microcontroller,  interface  and  a  laptop  computer.  The 
in.strument  guides  the  operator  in  the  collection  of  simultaneous  magnetic  field  and  position 
data  in  a  horizontal  plane  above  an  object.  Location  is  estimated  by  applying  a  custom- 
made  noniterative  location  estimation  algorithm  to  the  data.  A  byproduct  of  the  location 
algorithm  is  an  estimate  of  the  dipole  moment  which  is  used  by  a  continuous  parameter 
pattern  classifier  to  identify  the  object.  With  the  present  instrument,  6  to  13  minutes  arc 
required  to  collect  the  data,  location  and  dipole  moment  estimation  requires  5  seconds  and 
classification  requires  30  seconds. 

Detailed  experiments  were  perfonned  using  two  different  ferrous  spheroids,  similar  in 
shape  to  unexploded  ordnance,  in  order  to  determine  limits  of  error  in  estimation  of  the 
location  and  the  dipole  moment  and  the  error  in  classification. 

The  precision  in  estimating  location  varied  from  0.12  to  1.35%  of  the  depth  for  a  fixed 
object  and  orientation  at  a  given  depth.  There  was  a  slight  variation  in  location  estimate  with 
object  type  and/or  orientation  which  increased  the  location  uncertainty  to  1.67  to  1.90%  of 
the  depth.  The  precision  in  estimating  the  dipole  moment  varied  from  0.51  to  8.21  %  of  the 
dipole  moment  magnitude.  There  was  a  slight  increase  in  uncertainty  as  depth  decreased, 
but  the  precision  was  not  closely  correlated  with  peak  magnetic  field. 

All  three  components  of  the  location  estimate  had  a  depth  dependent  bias  which  was 
larger  than  can  be  accounted  for  by  the  precision  of  the  estimation.  The  Xq\  component 
of  the  location  was  almost  always  overestimated.  The  Xoz  component  was  generally 
overestimated  for  depths  less  than  0.55  m  and  underestimated  for  depths  between  0.55 
m  and  0.8m.  A  similar  trend  appeared  to  occur  for  the  X03  component  except  that  the 
crossover  from  overcstimation  to  underestimation  occured  at  about  0.45  m  depth.  Errors  in 
Xoi  were  generally  larger  than  those  for  Xoi  and  X(yi,  the  root  mean  square  (rms)  error  over 
all  depths  being  C,.)20  m  for  Xq\,  0.019  m  for  Xai  and  0.045  m  for  X03.  By  comparison, 
the  average  precision  for  a  fixed  object  and  orientation  was  0.005  m  for  Xoi,  0.003  m  for 
X02,  and  0.004  in  for  Xi^i. 

There  was  a  deviation  between  theoretical  values  and  experimental  estimates  of  the 
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dipole  moment  that  was  greater  than  the  uncertainties  in  the  estimates.  The  unweighted 
average  difference  between  estimated  and  theoretical  dipole  moment  components  as  a 
percentage  of  the  theoretical  dipole  moment  was  4.8  ±  7.6%  for  M\,  -3.3  ±  15.5%  for 
M2,  —10.8  ±  17.5%  for  and  24.5  ±  1 1 .4%  for  the  magnitude  of  the  vector  difference. 
There  was  no  clear  trend  with  object  type,  orientation  or  depth. 

Following  quantification  of  the  error  in  the  estimation  of  location  and  dipole  moment, 
sources  of  error  and  their  relative  effects  were  analysed.  It  was  argued  that  a  substantial 
fraction  of  the  uncertainty  (precision)  in  dipole  location  and  moment  vector  estimates  was 
due  to  positional  uncertainty  in  the  experiments. 

By  applying  the  location  algorithm  to  magnetic  total  field  data  generated  by  computer 
using  a  mathematical  model  based  on  magnetization  induced  in  spheroids  by  a  uniform 
magnetostatic  field,  the  effects  of  remnant  magnetization  and  positional  uncertainty  could 
be  eliminated.  It  was  shown  that  the  estimation  error  of  the  location  vector  was  roughly 
independent  of  whether  the  magnetic  field  data  was  experimental  or  computer  generated. 
Dipole  moment  estimation  errors  for  computer  generated  data  were  slightly  smaller  than 
those  for  experimental  data.  This  suggested  that  positional  error  in  these  experiments  had 
little  influence  on  location  and  moment  vector  estimation  and  that  remnant  magnetization 
in  the  spheroids  must  be  small  compared  to  the  induced  magnetization.  It  was  shown  that 
the  finite  volume  of  the  sensor  head  was  not  a  significant  source  of  error  for  the  geometry 
of  these  experiments.  A  variation  of  error  with  depth  and  object  orientation  was  noticed 
for  the  computer  generated  data.  This  suggested  that  the  main  contributor  to  the  error  was 
estimation  of  the  position  of  the  field  extrema  and  the  use  of  total  field  in  place  of  vertical 
component  magnetic  data. 

The  location  algorithm  was  applied  to  magnetic  vertical  component  field  data  generated 
from  the  computer  model  in  order  to  determine  how  much  of  the  error  in  estimation  was  due 
to  using  total  field  measurements  to  approximate  the  vertical  component  of  the  magnetic 
field.  Since  the  bias  in  estimation  of  the  sphere  for  vertical  component  data  was  roughly 
2  to  6%  of  the  moment  magnitude,  whereas  for  total  field  data  it  was  20%,  this  suggested 
that  the  error  due  to  the  total  field  approximati.m  varied  from  14  to  18%  of  the  moment 
magnitude.  It  was  estimated  that  the  error  due  to  e,stimation  of  extrema  positions  and  higher 
order  moments  varied  from  0  to  8%  of  the  moment  magnitude. 


The  analysis  of  error  in  the  computer  generated  total  magnetic  field  data  allowed  limits 
to  be  set  on  the  misclassification  rate  that  can  ultimately  be  attained.  Tlie  moment  estimates 
for  a  given  object  and  orientation  were  found  to  vary  by  between  2  and  4%  of  the  moment 
magrutude  as  the  depth  varied.  If  the  estimated  moments  were  used  in  the  design  set,  this 
suggests  that  the  limiting  misclassification  error  would  be  between  roughly  1.4  and  2.4% 
if  the  experimental  (mainly  positional)  error  were  substantially  less  than  the  algorithmic 
and  approximation  errors.  For  the  present  experimental  precision  (4.1  ±  3.0%  ),  the 
misclassification  error  should  be  roughly  between  2  and  5%  (assuming  quadrature  error 
summatioji)  if  the  design  set  were  experimentally  obtained. 

Pattern  classification  was  performed  with  a  computer  generated  dipole  moment  design 

UNCLASSIFIED  DRES-SR-582 


UNCLASSIFIED 


87 


set  consisting  of  8  objects  including  the  2  used  in  these  experiments.  The  gross  misclas- 
sification  rate  was  10.3%  for  the  class  number  4  spheroid  and  100%  for  class  number  6. 
Closer  examination  revealed  that  class  6  was  always  classified  as  class  number  3,  which  is  a 
similarly  shaped  spheroid.  This  is  likely  due  to  the  deviation  between  the  moment  estimates 
for  the  experimental  spheroid  and  the  computer  model.  The  repeatability  of  classification, 
that  is,  the  percentage  of  cases  in  which  an  object  is  classified  as  the  same  class,  is  a  better 
measure  of  classifier  performance.  TTie  repeatability  of  classification  was  91.7%.  If  we 
consider  only  the  two  objects  for  which  experimental  measurements  were  made,  we  see 
that  only  4/36  cases  (11.1%)  were  classified  incorrectly. 

The  overall  performance  of  the  smart  magnetometer  is  very  encouraging,  but  additional 
studies  need  to  be  performed.  Although  remnant  magnetization  has  been  shown  to  be 
small  in  these  and  previous  experiments,  this  must  be  verified  by  estimating  moments 
for  sets  of  identically  shaped  spheroids  at  fixed  depths.  This  should  also  be  repeated  for 
unexploded  ordnance.  The  possibility  of  accurately  determining  the  relative  contributions 
of  remnant  and  induced  magnetization  should  be  investigated.  A  major  source  of  error  is 
the  estimation  of  position  of  the  magnetic  field  extrema.  A  possible  improvement  is  to 
regularize  the  sampling  grid  for  the  magnetic  field  measurements  so  that  two  dimensional 
interpolation  may  be  used  to  estimate  the  extremum  position,  rather  than  its  centroid  as  is 
currently  done.  This  appproach,  which  will  be  investigated,  may  improve  the  accuracy  of 
the  peak  position  estimates.  The  possibility  of  applying  corrections  to  the  total  field  data  to 
improve  its  approximation  to  vertical  component  data  should  be  investigated.  One  possible 
method  to  achieve  this  might  be  to  use  a  three  sensor  orthogonal  axis  fluxgate  total  field 
magnetometer  which  is  either  on  a  stabilized  platform  or  whose  orientation  angles  with 
respect  to  a  fixed  coordinate  system  are  known.  Whether  such  a  system  is  practical  must 
be  investigated.  Most  importantly,  a  comprehensive  design  set  must  be  assembled  from 
measured  experimental  magnetic  field  maps  of  objects  of  interest. 

The  last  proposed  study  will  require  considerable  time,  given  that  between  6  and  13 
minutes  arc  required  to  collect  the  data  for  one  dipole  moment  estimation.  This  laboratory, 
in  conjunction  with  Pylon  Electronics  Ltd.  of  Ottawa,  Canada,  is  presently  developing  a 
multisensor  fluxgate  gradionieter  which  will  incorporate  a  modified  Wynn-Frahm  location 
algorithm  [16]  and  the  present  pattern  classifier.  Ihe  gradiometer  will  be  able  to  estimate 
location  from  a  single  measurement  in  space  and  hence  sufficient  data  for  a  moment 
estimation  will  be  collected  in  less  than  a  second.  To  use  a  design  set  from  the  gradiometer 
in  the  magnetometer,  the  estimated  dipole  moments  from  one  instrument  must  be  nearly 
equal  to  those  from  the  other.  A  comparative  study  of  the  dipole  moment  estimation  errors 
for  both  the  smart  magnetometer  and  gradiometer  will  be  conducted  using  the  spheroid 
model,  and  later  using  real  magnetic  field  and  field  gradient  data. 
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